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Abstract 

In this review, analytical results obtained for a wide class 
of stationary axisymmetric flows in the vicinity of com- 
pact astrophysical objects are analyzed, with an empha- 
sis on quantitative predictions for specific sources. Re- 
cent years have witnessed a great increase in understand- 
ing the formation and properties of astrophysical jets. 
This is due not only to new observations but also to 
advances in analytical theory which has produced fairly 
simple relations, and to what can undoubtedly be called 
a breakthrough in numerical simulation which has en- 
abled confirmation of theoretical predictions. Of course, 
we are still very far from fully understanding the phys- 
ical processes occurring in compact sources. Neverthe- 
less, the progress made raises hopes for near-future test 
observations that can give insight into the physical pro- 
cesses occurring in active astrophysical objects. 
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1 Introduction 

This review was not specially written for the Uspekhi 
Fizicheskikh Nauk (Physics- Uspekhi) issue devoted to 
the memory of V L Ginzburg. Nevertheless, I would 
like to hope that the spirit of this review is close to that 
of other papers in this issue. I was happy to work closely 
with Vitaly Lazarevich for more than 30 years, starting 
from my graduate student days, and he definitely played 
a significant role in my scientific development. So, this 
review, hopefully, bears a fraction of the soul of Vitaly 
Lazarevich. 

Vitaly Lazarevich was a passionate person. Astro- 
physics, undoubtedly, fascinated him most. But the 
scale of his personality was such that this passion did 
not separate, but instead united, people. So, it is not 
surprising that at the Lebedev Institute the astrophysi- 
cal seminar headed by Vitaly Lazarevich for more than 
several decades continues, and scientists from many in- 
stitutes participate in its work. The Department of 
Physics and Astrophysics Problems at Moscow Institute 
of Physics and Technology (MIPT), which Vitaly Lazare- 
vich founded in 1968 and headed until recently, continues 
to be one of the leading institutions in teaching young 
astrophysicists. Numerous pupils of V L Ginzburg and 
pupils of his pupils working in the leading astrophysical 
centers of the world keep his unique trademark in their 
studies. 

The astrophysical heritage of V L Ginzburg is enor- 
mous. He obtained fundamental results in the theory of 
propagation of electromagnetic waves in cosmic plasma. 
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in the theory of the origin of cosmic rays, and in the 
theory of neutron stars and black holes. In all cases, a 
simple model allowing the understanding of the essence 
of physical process in observed astrophysical sources laid 
the basis of the theory. The present review, hopefully, 
was written in the same spirit. 

Astronomy, as follows from the very appellation, is 
the science that stemmed from the observations of stars. 
During hundreds of years the people observed stars in the 
sky and gained insight into the laws of Nature. The stars 
appeared to be always unchanged and existing for an in- 
finite amount of time. After the appearance of spectral 
analysis, the first astrophysical observations and, later, 
the theory of radiation generally confirmed this point of 
view. The lifetimes of most stars turned out to be com- 
parable to the age of the Universe. Thus, in the 1950s, 
when radio astronomy began, stars emitting thermal ra- 
diation seemed to be the main objects for studies. In ra- 
dio astronomy, the brightness temperature remains even 
now the basic characteristic of radiation intensity. 

However, the first radio astronomical observations, 
and especially observations in the X-ray and gamma-ray 
ranges, which started in the middle of the 1970s, dis- 
covered numerous nonthermal sources in the Universe. 
These objects are sufficiently compact (i.e., the spatial 
resolution of the existing detectors is insufficient to de- 
termine their internal structure) and, in addition, are 
highly variable. In active galactic nuclei the variabil- 
ity timescale (months or sometimes even days) is small 
according to the cosmic timescale, with the variability 
timescale of radio pulsars and sources of gamma-ray 
bursts being the fractions of a second, which is small 
even to Earth's measures. The activity, i.e., high vari- 
ability on timescales r ~ R/c, as well as the generation 
of nonthermal radiation indicates that in most cases we 
are dealing with relativistic objects, namely, with objects 
in which matter moves with velocities close to the speed 
of light. 

Jet eruptions represent one of the visible appearances 
of the activity of compact astrophysical objects. We shall 
briefly discuss their properties in Section 2. They arc ob- 
served in both relativistic objects (such as active galac- 
tic nuclei and microquasars) and in young stars where 
the motion of matter is definitely nonrelativistic. This 
means that we are dealing with some universal and ex- 
tremely efficient mechanism of energy release. There- 
fore, the key theoretical problems include the question of 
the energy source of the activity of compact objects, the 
understanding of their energy release mechanism, and 
the collimation of matter outflows. We shall postpone 
the detailed discussion of arguments against alternative 
models until the next section, and here we only remind 
the main arguments favoring the magnetohydrody-namic 
model of activity of compact sources, which is accepted 
by most astrophysicists. 

The model of the unipolar inductor, i.e., the source 
of direct current, lies at the heart of the magnetohy- 
drodynamic approach. As we shall show in Section 2, 
conditions for the existence of such a 'central engine' 
are satisfied in all the compact sources discussed below. 



Indeed, all compact sources are assumed to harbor a 
rapidly spinning central body (black hole, neutron star, 
or young star) and some regular magnetic field, which 
leads to the emergence of strong induction electric fields. 
The electric fields, in turn, lead to the appearance of 
longitudinal electric currents and effective particle ac- 
celeration. The collimation mechanism in this model is 
related to the well-known property of mutual attraction 
of parallel currents. 

The first studies of the electromagnetic model of com- 
pact sources (namely, radio pulsars) were carried out as 
early as the end of the 1960s [1-4]. It was evidenced that 
there are objects in the Universe in which electrodynam- 
ical processes can play the decisive role in the energy re- 
lease. Then, in 1976 R Blandford [5] and R Lovelace [6] 
independently suggested that the same mechanism can 
also operate in active galactic nuclei. In the same year, 
G S Bisnovatyi-Kogan, Yu P Popov, and A A Samokhin 
proposed a magnetorotational mechanism of the super- 
nova explosion [7] (i.e., the model of an essentially non- 
stationary phenomenon), in which jet eruptions can also 
be formed [8] . This model has remained the leading one 
for nearly 40 years. However, only recently have some 
key properties become clear. This is related both to 
advances in the theory which have at last formulated 
sufficiently simple; analytical relations, and to the break- 
through in numerical simulations which confirmed theo- 
retical predictions. 

The reader can find the detailed introduction to the 
analytical theory in the author's monograph [9] (see also 
the review in Physics-Uspekhi in 1997 [10]). However, 
first, the monograph was devoted to the basics of the 
theory, and qualitative predictions for specific astrophys- 
ical sources were discussed only very briefly. Second, 
the monograph clearly could not include the results of 
numerical calculations carried out in the last five years 
since its publication. This is the main reason for writing 
the present review. In addition, here we shall correct 
formulas from the monograph in which misprints were 
found. 

Of course, we are still far away from the full under- 
standing of the essence of physical processes proceeding 
in compact sources. In fact, now we have only agreement 
between theory and numerical modeling. All results have 
been obtained applying ideal one-liquid magnetohydro- 
dynamics, though by different methods (the theory is 
based on stationary equations, while numerically the 
time relaxation problem is solved). In particular, it is not 
yet clear which of the main physical characteristics of the 
central engine (such as the mass of the central body or 
its rotation velocity) should fully determine the observed 
energy release. Nevertheless, the progress achieved over 
recent years raises hopes for test observations already in 
the nearest future, which can give insight into physical 
processes occurring in active astrophysical sources. 



MHD MODELS OF ASTROPHYSICAL JETS 



3 



2 Jets 

2.1 Active galactic nuclei 

The main properties of the central engine in active 
galactic nuclei, which are presently accepted by most 
astrophysicists, can be summarized as follows [11,12]. 
In the center of the host galaxy there is a supermas- 
sive black hole (its mass reaches 10^-10^ M0, where 
Mq « 2 X 10'^'^ g is the mass of the Sun), onto which 
accretion of the surrounding matter occurs [13]. Only in 
this case it is possible to explain the very high efficiency 
of the energy release and the compactness of the central 
engine. The energy source of activity of galactic nuclei 
can be related to both the rotational energy of the black 
hole, viz. 

and the energy of the accreting matter. Here 
Tg = 2GM/c^ is the radius of the black hole, Jr is 
the moment of inertia, fin and M are the angular ve- 
locity and the mass of the black hole, respectively, and 
c is the speed of light. The existence of supermassive 
objects is also supported by the fact that the Eddington 
luminosity 

(i.e., the luminosity at which the gravitational force act- 
ing on the accreting matter is balanced by the radiation 
pressure force) is close to the characteristic luminosity of 
active galactic nuclei [14]. Moreover, the duration of the 
active phase td = -BtotZ-^^Edd estimated using formulas 
(1) and (2) is on the order of 10^ years, which is also in 
agreement with observations. 

Further, it is usually assumed that the accretion of 
matter proceeds through a disc [15]. Thus, the prefer- 
ential direction — the axis of rotation — emerges natu- 
rally in space, along which the formation of jets is possi- 
ble. As a black hole itself cannot have the self-magnetic 
field (the so-called 'no-hair theorem'), the generation of 
a large-scale magnetic field in the vicinity of the black 
hole is believed to occur in the accretion disc [16-18]. 

According to the modern concept, massive central ob- 
jects are present in most galaxies and remain active 
only if a sufficient amount of matter falls on them. 
This restricts their active lifetime. Unfortunately, as 
stated above, the angular resolution of modern de- 
tectors does not allow us to directly observe plasma 
flow on the scales comparable to the black hole size 
rg « 3 X 1O"(M/1O^M0) cm. Therefore, we have to 
judge the activity of galactic nuclei only using indirect 
evidence, by observing flows on much larger scales. 

Let us remember that the diffuse radio emission 
around active galaxies is observed from regions located 
at distances of tens or even a hundred kiloparsecs from 
their nuclei. Very shortly after the discovery of these 




Figure 1: Radio image (5 GHz) of active regions and 
jet eruptions from the nucleus of the Cygnus A galaxy 
[20]. The distance between bright spots is about 80 kpc, 
which is 9-10 orders of magnitude greater than the size 
of the central black hole. 
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Figure 2: Radio image of the jet eruption from the galaxy 
M87 near the central engine [22]. The jet transverse size 
is about 1 pc. 



regions at the beginning of the 1960s, this emission was 
associated with coUimated plasma ejections (jets) flow- 
ing out the galactic nuclei [12]. It is precisely these jets 
that transport matter and energy from the active nuclei 
to those regions (Fig. 1). Observations show that the 
jets can be accelerated and collimated very close to a 
galactic nucleus. For example, in the case of the near- 
est active galaxy M87 the formation of the jet occurs 
within a radius of 60 Tg from the nucleus [19]. In recent 
years, the internal structure close to the jet base was 
resolved in several sources, where the jet tranverse di- 
mension usually does not exceed several parsecs [21, 22] 
(Fig. 2). ^ ^ 

The matter in jets from active galactic nuclei has a 
very high energy — the bulk Lorentz factor of a jet is at 
least a few unities. For example, this motion is directly 
observed in the M87 galaxy, with the bulk Lorentz fac- 
tor of the outflow being 7^6 [23]. In many cases, 
the matter continues moving with relativistic velocities 
up to huge distances from the nucleus before noticeable 
braking due to interaction with the ambient intergalactic 
medium. Another peculiar feature of jets is their high 
degree of coUimation within a cone characterized by an 
opening angle of only several degrees. 

Unfortunately, observations do not yet allow reliable 
estimations of the energy and mass fluxes in jets from 
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active galactic nuclei, of the magnitude of the magnetic 
field both close to the black hole and in the jet itself, or of 
the composition of jet eruptions. The spectrum of radi- 
ation from galactic jets (in contrast, for example, to the 
spectrum of jets from young stars) does not exhibit any 
spectral features of moving matter, i.e., neither atomic 
(ionic) lines nor the electron-positron pair annihilation 
line are observed. To this regard, there are arguments 
both in favor [24] and against [25] the leading role of 
electron-positron plasma, so it is now impossible to say 
exactly which mechanism of energy transfer to the jet 
actually operates. 

Where a physical nature of the galactic nucleous ac- 
tivity is concerned, several mechanisms of particle ac- 
celeration and jet collimation have been proposed, but 
so far there is no definite answer as to which of them 
are actually realized. It is possible that different mecha- 
nisms operate in different sources, or, just the opposite, 
all mechanisms are realized simultaneously. 

Gas- dynamic acceleration. The acceleration and colli- 
mation of a jet can be related to the presence of an ambi- 
ent medium with high pressure which decreases with dis- 
tance from the center [26, 27]. Such a medium could play 
the role of an external wall collimating the outflow. The 
pressure of the external hot medium can, in principle, 
be estimated from X-ray observations [28]. This mecha- 
nism possibly explains how weak jets in our Galaxy and 
in some Seyfert galaxies (i.e., low- active galaxies) are 
formed. On the other hand, the observed pressure of the 
hot matter aroimd the most powerful jets from active 
galactic nuclcii is not siifficiently high, and there miist be 
an alternative mechanism of plasma confinement. 

Acceleration by radiation. As the photon density near 
the central source can be very high, the radiation-driven 
mechanism of jet matter acceleration by radiation pres- 
sure was proposed [29, 30]. In this model, it is assumed 
that the inner parts of the disc can serve as a nozzle 
directing matter outflows accelerated by the radiation 
pressure. However, this mechanism also meets some dif- 
flculties. For example, there is no correlation between 
the jet power and the luminosity of the source — many 
sources with very powerful jets are low-luminous sources 
[31]. Another diflSculty comes from the fact that, start- 
ing from the sufliciently low particle energies 7 ~ 3, the 
radiation field more effectively brakes particles than ac- 
celerates them [32]. This contradicts observations of 'su- 
perluminal' sources in which the energy of plasma par- 
ticles is much higher. In addition, if the jet was formed 
in a system with a thin accretion disc emitting radiation 
more or less isotropically, additional mechanisms for the 
jet collimation should be invoked. A modification of this 
model involving the formation of a funnel in a thick ac- 
cretion disc can explain the initial jet collimation, but 
there are indications that such a structure is unstable 
[31]. 

Magnetohydrodynamic mechanism. As noted above, 
most researchers favor the magnetohydrodynamic model 
of jet formation. The magnetohydrodynamic (MHD) 
model was successfully utilized to describe many pro- 
cesses in active nuclei, and, in particular, in connection 



with the problem of the origin and stability of jets, as 
well as to explain the energetics of processes proceeding 
near the central black hole. The magnetic field here is 
the natural link between the central engine and the jet. 
Moreover, in this model it is easy to understand why the 
jet matter can predominantly consist of electron-positron 
plasma. As was shown in Refs [33-35], it can be gener- 
ated on the magnetic field lines threading the black hole 
horizon. 

In the simplest version, the picture is as follows: the 
regular magnetic field generated in the disc links the 
spinning central engine (the disc and the black hole) 
with infinity. The plasma outfiow occurs along the mag- 
netic field lines; the electromagnetic energy flux is also 
directed along the magnetic field lines. The longitudi- 
nal electric current flowing along the jet forms a toroidal 
magnetic field, and the magnetic field pressure associ- 
ated with this toroidal component can coUimate the jet. 

It should be noted, however, that in a real astrophys- 
ical system the total current flowing from the central 
engine should vanish, so the Ampere force in the current 
closure region will, on the contrary, decollimate the flow 
(the antiparallel currents repulse). Therefore, an exter- 
nal medium (for example, a subrelativistic wind outflow 
from the accretion disc) is necessary to coUimate the jet. 
In addition, the question as to whether it is possible to 
consider a black hole immersed in the external magnetic 
field as a unipolar inductor turned out to be also rather 
nontrivial. It required almost 30 years of studies after 
the paper by Blandford and Znajek [33], which laid the 
basis of the theory in 1977, before the needed clarity was 
reached in this question. We shall discuss these points 
in more detail in Section 3. 

2.2 Microquasars 

Microquasars comprise galactic objects in which the jet 
formation is due to accretion onto a compact relativis- 
tic object (neutron star or black hole). In other words, 
all microquasars reside in sufficiently close binary sys- 
tems in which the effective flow of matter from the star 
companion occurs. The rate of matter inflow in such 
systems is larger than can be swallowed by the central 
object. As a result, some accreting matter that carries, 
in particular, an excessive angular momentum is expelled 
from the system in the form of jets. Observations of mi- 
croquasars show that jets are related to thick accretion 
discs. In other words, no jets are known for systems with 
thin discs. The reason for that is unclear: either a thin 
disc insufficiently coUimates the outfiow, or the magnetic 
field generated by the thin disc is not strong enough. 

Microquasars represent a small population of ob- 
jects, including only around ten sources [36], with only 
half of them demonstrating noticeable relativistic jets 
{v > 0.9 c). The characteristic longitudinal size of jets is 
usually 0.1 pc, with the jet spread angle being within 
several degrees (Fig. 3). The total energetics are about 
]^q37 gj-g g-i ^j^g relativistic velocity of the bulk 

motion of matter in the jets, some sources demonstrate 
the super luminal motion effect, with the apparent an- 
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Figure 3: Radio image of the jet from the microquasar 
iE1740. 7-2942. The characteristic length of the jet mea- 
sures 0.1 pc [38]. 



gular velocity being several orders of magnitude larger 
(due to a relative proximity of these objects) than that 
observed in jets from active galactic nuclei. 

Historically, the first revealed object of this class was 
the famous source SS433 [37] in which, however, the gas 
ejection velocity in jets is only 0.26 c. Such a velocity 
can be easily explained by the radiation pressure from 
highly heated internal regions of the accretion disc. As 
for relativistic jets, the first source was discovered only 
in 1994 [38]. Since the appearance of near-light velocities 
due to radiation or gas pressure is problematic, it has not 
been ruled out that, to explain them, an electrodynamic 
model similar to that used in explaining the origin and 
coUimation of extragalactic jets should be invoked again. 
This model is also supported by the fact that in all but 
one microquasar (SS433) no emission lines from jets are 
observed. This indirectly points to the electron-positron 
composition of matter in jets [36]. Finally, it should 
be noted that in most microquasars the jet is separated 
in individual blobs at large distances from the central 
engine, which is thought to be due to a long duty cycle 
of the work of the central engine. 

2.3 Sources of cosmological gamma-ray 
bursts 

As regards the sources of cosmological gamma-ray 
bursts, there are indirect, although sufficiently reliable, 
arguments in favor of the presence of jets related exactly 
to relativistic strongly magnetized outflows, which we 
shall discuss in this review. It is well known that 
the discovery of the optical afterglow [39], as well as 
afterglows in other spectral ranges, which allowed the 
measurement of the distance to these sources from 
the observed redshifts of the host galaxies, put serious 
constraints on their energetics [40]. If the observed 
gamma-ray radiation were emitted isotropically, the 
total energy release for the typical distance to these 
sources of several gigaparsecs would reach 10^'' erg. 
However, we do not know at present processes with 
such huge energy liberation. On the other hand, the 
small duration of the burst 10 s) restricts the size of 
the emission region, which, in turn, does not allow us 
to explain the observed nonthermal gamma-ray spectra. 



since the optical depth in the source proves to be very 
high [41]. 

If it is assumed that gamma-rays are emitted within 
a narrow cone angle -d 1°, the observed energy can 
be reduced to 10^^ erg, which is already to an order of 
magnitude of the energy release during supernova explo- 
sions. On the other hand, the observed optically thin 
nonthermal gamma-ray spectra immediately imply the 
presence of ultrarelativistic outflows with bulk Lorentz 
factors of ^ 100-300. Only in this case can the com- 
pactness problem of the source be resolved, since the 
estimated size of the emitting region also increases as 
the square of the bulk Lorentz factor (i.e., by 10'*-10^ 
times), and the optical depth, which is proportional to 
the density multiplied by the size of the region, decreases 
respectively by 10*-10^° times. 

However, the ultrarelativistic character of the outflow, 
in turn, puts constraints on the particle composition in 
the jet, since the presence of a significant fraction of 
baryons with such energy in the outfiow would contra- 
dict the total energy release in the gamma-ray burst. 
Therefore, the contribution of protons must be smaller 
than 10~^ of the total number of particles, so that only 
electron-positron jets should be considered. The exis- 
tence of jets is also evidenced by the presence of the 
characteristic bend of the light curve of the afterglow, 
when the power law index a in the radiation intensity 
dependence on time, Wtot oc t~", changes from a ~ 1.1 
to a ~ 2.0 after a span of about a few days following the 
burst. This effect is related to the cessation of relativistic 
contraction of the radiation cone in the direction of parti- 
cles' motion toward an observer. Incidentally, this model 
allowed independent confirmation of the jet spread angle 
z9 - 1° and the bulk Lorentz factor 7 ~ 100-300 [42]. 

The nature of the central engine giving rise to strongly 
magnetized jets can be usually related to the collision of 
two neutron stars [43, 44] or of a neutron star and a 
black hole [45], or, most likely, to the collapse of the 
massive core of an unusual supernova [46, 47]. How- 
ever, in most models a rapidly spinning solar-mass black 
hole ultimately emerges, which loses its rotation energy 
via the Blandford-Znajek process [45, 48 50]. Indeed, as 
we have seen, this process easily provides a natural ex- 
planation for both the low baryonic load of the jet and 
the large bulk Lorentz factors of jet particles. In other 
words, the model again is constructed similarly to the 
scheme proposed for active galactic nuclei. In partic- 
ular, the key processes here also include the magnetic 
field generation in the plasma around the black hole, the 
interaction of the black hole with the accretion disc via 
magnetic field lines, and the generation of particles in the 
magnetosphere. To explain the observed energy release, 
it is necessary to assume that the magnetic field near 
the black hole must be as high as lO^** or even 10^^ G. 
The generation of such a high field is thought to be pos- 
sible in nonstationary processes like the supernova core 
collapse or binary neutron star coalescence [51, 52]. 
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2.4 Radio pulsars 

The discovery of radio pulsars at the end of the 1960s, 
which axe the sources of pulsating cosmic radio emis- 
sion with the characteristic period P ~ 1 s [53], is def- 
initely one of the major astrophysieal discoveries of the 
20th century. Indeed, for the first time a new class 
of cosmic sources related to neutron stars, whose ex- 
istence was theoretically predicted away back in the 
1930s [54], was discovered. Neutron stars (mass of about 
1.2-1.4 Mq, and radius R of only 10-15 km) must result 
from the catastrophic compression (collapse) of usual 
massive stars at the late stage of their evolution or, for 
example, of white dwarfs whose mass exceeds the Chan- 
drasekhar mass limit of 1.4 Mq due to accretion from 
the companion star. It is this formation mechanism that 
provides the simplest explanation for both small spin pe- 
riods P (the smallest known spin period P = 1.39 ms) 
and superstrong magnetic fields with Bq ~ 10^^ G [1,2]. 

Interestingly, the basic physical processes determining 
the observed activity of radio pulsars were understood 
almost immediately after their discovery. For example, 
it became clear that highly regiilar pulsations of observed 
radio emission are related to the rotation of neutron 
stars. Next, radio pulsars are powered by the rotational 
energy of the neutron star, and the mechanism of energy 
release is related to their superstrong magnetic field with 
Bo ~ 10"'^^ G. Indeed, energy losses estimated using the 
simple magnetodipole formula [44] are as follows: 

W,,, = -jM^l^^^sm''x, (3) 
6 cr 

where Jp ^ Mi?^ is the moment of inertia of the neu- 
tron star, X is the angle between the magnetic dipole 
axis and the spin axis, and fl = 27r/P is the angular 
velocity of the neutron star rotation. For most pulsars, 
energy losses range from 10"^^-10"^** erg s~^. These energy 
losses exactly correspond to the observed spin-down rate 
dP/dt ^ 10""'^^, or to the spin-down time td = P/P ^ 
1-10 mln years. Let us keep in mind that the frac- 
tion of radio emission amounts to only 10^^ 10^® of 
total energy losses. For most pulsars this corresponds 
to 10^^-10^^ erg s~^, which is 5-7 orders of magnitude 
less than the luminosity of the Sun. 

As shown in Refs [56, 57], the actual energy losses 
cannot be due to magnetodipole radiation because the 
plasma that fills the magnetosphere will fully screen the 
low-frequency radiation from the neutron star. However, 
energy losses can be caused by longitudinal electric cur- 
rents circulating in the magnetosphere and looped across 
the surface of the central engine. As a result, in this case, 
too, the main energy release near the neutron star is re- 
lated to the electromagnetic energy flux (the Poynting 
vector flux), and the total energy losses can be again 
estimated using formula (3). 

Most radio pulsars constitute single neutron stars. Of 
the 1880 pulsars known by the middle of 2010, only 140 
were members of binary systems. However, in all these 
cases it is reliably known that there is no somewhat ap- 
preciable mass transfer from the companion star to the 



neutron star. Since, as already stressed, the radio lu- 
minosities of pulsars are low, the modern sensitivity of 
detectors allows observations of pulsars only up to dis- 
tances of 3 5 kpc, which is smaller than the distance to 
the galactic center. Therefore, we can observe only a 
small fraction of all 'active' pulsars. The total number 
of neutron stars in our Galaxy must be around 10^-10^. 
Such a big number of extinguished neutron stars can be 
naturally related to the small duration of their active 
life, as discussed above. 

The jets are only observed in Crab and Vela radio 
pulsars [58, 59], which is not surprising, since, in con- 
trast to the compact objects considered above, the pulsar 
magnetosphere is not axisymmetric. On the other hand, 
only axisymmetric configurations were actually consid- 
ered until recently in the theory of pulsar wind. Based 
on these studies, the main features of strongly magne- 
tized winds were understood. Nevertheless, even in this 
approximation for smooth flows, it has thus far been 
impossible to construct a self-consistent model which 
jointly describes the energy transfer from the neutron 
star surface to infinity and includes effective particle ac- 
celeration, i.e., an almost complete transformation of 
the electromagnetic field energy into the energy of the 
plasma flowing out. Because of this, different models 
are actively being discussed at present, which, to vari- 
ous degrees, propose going beyond the framework of the 
'classical' scheme (see, for example, Refs [60-62]). 

Indeed, observations show that most energy far from 
the neutron star must be carried by relativistic parti- 
cles. For example, the analysis of the emission from the 
Crab Nebula in the shock region located at a distance of 
~ 10^^ cm from the pulsar in the region of interaction of 
the pulsar wind with the siipernova remnant definitely 
shows that the total flux Wcm of the electromagnetic en- 
ergy in this region is no more than ^ 10^^ of the particle 
energy flux Wpart [63]. Thus, the Poynting vector flux 
in the asymptotically remote region must be completely 
converted into the outgoing plasma flux. The presently 
known axisymmetric numerical models of jets from ra- 
dio pulsars [64-66] were constructed exactly under this 
assumption. 

However, the transformation apparently occurs al- 
ready much closer to the neutron star, namely at dis- 
tances comparable to the size of the light cylinder. This 
is evidenced by the detection of variable optical emission 
from companions in some close binary systems involving 
radio pulsars [67]. This variable optical emission with a 
period equal exactly to the orbital period of the binary 
can be naturally related to the heating of the compan- 
ion's part facing the radio pulsar. It was found that the 
energy reradiated by the companion star almost matches 
the total energy emitted by the radio pulsar into the 
corresponding solid angle. Clearly, this fact cannot be 
understood either in the magnetodipole radiation model 
or by assuming a Poynting-dominated strongly magne- 
tized outflow, since the transformation coefficient of a 
low-frequency electromagnetic wave cannot be close to 
unity. Only if a significant fraction of the energy is re- 
lated to the relativistic particle flux can the heating of 
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Figure 4: Optical image of jets from the system HH47 
(see, for example, Ref. [70]). The scale corresponds to 
1000 a.u. 



the star's surface be effective enough. Therefore, the so- 
called <T-problem — the question as to how the energy is 

transferred from the electromagnetic field to particles in 
the pulsar wind — remains one of big puzzles in modern 
astrophysics. 

2.5 Young stars 

Jets from young stars were indirectly discovered at the 
begirming of the 1950s, when G Herbig and G Haro [68, 
69] discovered a new class of extended diffuse objects 
usually existing in pairs and, as became clear later, con- 
nected by thin jets with young rapidly rotating stars [70] . 
The formation of such jets can naturally be related to the 
need of removing most effectively the excessive angular 
momentum that prevents the formation of a star. As 
we see, the situation here is quite similar to that with 
active galactic nuclei, where first a diversity of differ- 
ent types of sources (quasars, Seyfert galaxies, and radio 
galaxies) were discovered, and only later on did it be- 
come clear that the activity of all these sources has a 
similar nature. Moreover, the similarity of the observa- 
tional features suggests that the physical mechanism of 
jet formation from young stars can also be similar to that 
from active galactic nuclei. And this is despite the fact 
that physical conditions near a young star (mass of order 
3—10 Mq, and total energy release ranging from 10^^ to 
Iq36 j,j.g g-1 g^j.^, dramatically different from those in the 
centers of active galactic nuclei. One of the main differ- 
ences here is the nonrelativistic character of gas outflow 
from young stars. 

Presently, more than 250 Herbig-Haro objects are 
known [71]. As shown in Fig. 4, they represent bright 
condensations with an angular size of several seconds 
of arc (linear size of order 500-1000 a.u.), usually sur- 
rounded by a bright diffuse envelope. Their spectra 
mainly show emission lines of hydrogen and some other 
low-excitation elements. A shock wave propagating with 
velocities 40-200 km s~^ through a gas with a density of 
~ 10^ cm~^ is apparently the main source of excitation 
[70]. 

As in the case of radio galaxies, the activity of Herbig- 
Haro objects is dictated by collimated outflows which are 
well seen in forbidden lines. Nearly 60% of the objects 
demonstrate both jets, while in other cases the receding 
jet is blocked by the accretion disc. The extent of the op- 
tical jets is of order 0.01-2 pc, and their velocity reaches 




I 1 HH30 

Figure 5: Formation of a jet from a young star in the 
system HH 30 [70]. The accretion disc is clearly seen. 
Here also the scale corresponds to 1000 a.u. 



600 km s''^. The gas density in the jets is estimated 
to be 10-100 cm^'^, and the mass outflow rate comes 
to 10~^-10~^° Mq yr~^. The degree of collimation of 
the jets (the ratio of the observed length to the width) 
can be as high as 30. The total jet opening angle is in 
the range of 5-10°. In addition to highly elongated jets, 
molecular outflows with a much smaller collimation de- 
gree are observed near young stars. Their size may run to 
0.04 4 pc, and the velocity of gas motion does not exceed 
5-100 km s~^. Here we should stress that this velocity 
is much higher than the speed of sound in an outflow 
with a temperature of only 10- 90 K. The total mass of 
the ejected gas is estimated to be 0.1 200 Mq, and the 
total kinetic energy stored in the molecular outflows can 
reach 10^^ and even 1047 erg. The direct observation of 
rotation of the jets is the most important recent discov- 
ery. The characteristic velocities at an axial distance of 
20-30 a.u. range from 3-10 km s^^ [72, 73]. There is 
also direct evidence of the spiral structure of the mag- 
netic field in the jets [74] . All these facts unambiguously 
support the MHD model. 

As in jets from microquasars, a strong instability fre- 
quently develops in collimated outflows from young stars 
at large distances from the central engine (see Fig. 4), 
so that the outflow is split into separate blobs. On the 
other hand, as seen from Fig. 5, the flow near the base 
of the jet can be considered sufficiently regular. 

As for the physical nature of collimated jet formation, 
this question is still far from solved. It is only clear 
that the power of the central engine is always sufficient 
to accelerate the outflowing gas; however, the mech- 
anism of energy transformation remains unclear. We 
stress that in contrast to relativistic galactic objects (for 
example, microquasars), where the formation of jets is 
possibly caused by supercritical accretion, the luminos- 
ity in yomig stars never approaches the Eddington limit. 
On the other hand, it is clear that the key role in the 
collimated outflow formation is just played by accretion 
discs which undoubtedly exist around young stars. This 
is supported by the direct correlation between the power 
of the gas flux and the mass of the disc, estimated from 
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its luminosity, as well as some other correlations [75, 76]. 
The parameters of the discs can be very different. For 
example, their masses range from 0.1-100 Mq, while the 
outer radii can vary from 10 a.u. to 0.1 pc. 

It is important that, in contrast to discs around rel- 
ativistic objects (neutron stars and black holes), the 
gas temperature in discs around young stars is only 
20-100 K. As a result, as in the case with active galac- 
tic nuclei, neither the radiation pressure force nor gas 
pressure can explain the high velocities observed in the 
collimated outflows [71]. Therefore, to explain the jet 
formation and particle acceleration, models in which the 
magnetic field plays the key role and effectively medi- 
ates the interaction between the accretion disc and the 
jet were invoked once again. Because the real structure 
of the magnetic field in the proximity of a young star 
is presently unknown, here, too, both models in which 
the magnetic field of the star itself has a dominant role 
[77] and models in which the magnetic field of the disc 
plays the decisive role [75, 78] have been proposed. It is 
seen that here we meet the same problems regarding the 
structure of the initial magnetic field as in the study of 
the black hole magnetosphere. 

3 Basics of the MHD approach 

3.1 The key idea — unipolar inductor 

As already said, the notion of a unipolar inductor is the 
main physical idea that underlies the MHD theory of 
compact objects. Referring to Fig. 6, a rotating magne- 
tized ball can serve as the battery that determines the 
energy release from the central engine. Indeed, assuming 
the high conductivity of the ball, the freezing-in condi- 
tion of the magnetic field, viz. 

ri X r 

Ein + X Bin = (4) 

c 

(i.e., simply the condition that the electric field in the 
rotating reference frame vanishes), leads to the appear- 
ance of the potential difference SU between points a and 
b. To an order of magnitude, this potential difference 
can be 

SU ^ ERo ^ ^B, (5) 
c 

where Rq is the transverse size of the working area. As a 
result, the total energy release Wtot on the external load 
TZ will be given as 

Wtot = ISU, (6) 

where the electric current / = 6U/TZ. Here, however, 
several conditions should be met. First, the electric 
circuit must touch the ball at different latitudes, i.e., 
at points with different electric potentials. Second, the 
electric circuit should rotate with an angular velocity 
different from that of the magnetized ball. The current 
flowing along a wire tightly welded on the ball will be 
absent. 



Load 1^ 




Figure 6: The unipolar inductor as the source of a direct 
current. Inside the magnetized ball, the electric current 
flows against the electric field direction. 



We stress that the energy source [electromotive 
force (EMF)] in the unipolar inductor is due to the 
kinetic energy of rotation. Indeed, as seen from Fig. 
6, charges inside the ball move against the direction 
of the electric field. This becomes possible due to the 
force by which the lattice acts on charges carried along 
the wire, which violate the freezing-in condition inside 
the ball. Conversely, the Ampere force acting from the 
side of the surface electric current on the ball's material 
brakes its rotation. Therefore, the principle of work 
of the unipolar inductor (or, as it is sometimes called, 
the unipolar Faraday generator) is not the Faraday 
effect as such (where the EMF induced in a current 
loop depends on the variation of the magnetic flux), 
since the flux through the circuit remains constant. 
Notice that the reverse situation is also possible: if 
one applies a potential difference to a magnetized 
ball (i.e., if one replaces the load in Fig. 6 by the 
voltage source), the ball starts rotating. On the site 
http : / /fiziks . org^ua/sainyj-prosto j-y-inire-elektrodvigatel /[ 
which is devotlxf to laboratory studies in secondary 
school, one can find a video illustrating the work of such 
a device. 

As we have understood, for the central engine to op- 
erate it is necessary to have: 

• rotating body; 

• regular magnetic field, and 

• well-conducting wire. 

Then the current, and hence the energy losses, will be 
determined by the value of the external resistance TZ. 
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Let us see now whether these conditions can be met in 
compact astrophysical objects. 

As we have seen, a central rotating body in active 
astrophysical sources is undoubtedly present. For ex- 
ample, the spin periods of young stars are about sev- 
eral days (the inner parts of accretion discs rotate even 
faster). The spin periods of most radio pulsars are close 
to 1 s; however, they can be as small as a few mil- 
liseconds, which is already close to the limiting speed 
of rotation {flR/c ~ 0.1). The rotational velocities of 
black holes in active galactic nuclei, to tell the truth, are 
unknown, but we can suppose that due to disc accre- 
tion (it is in this way that the millisecond-period pul- 
sars are thought to have been spun up) their spin pa- 
rameter r^H-R/c = a/2M (see the Appendix) can also 
be sufficiently large. For example, the estimate of the 
black hole rotational velocity in the nucleus of Seyfert 
galaxy MCG 06-30-15, as inferred from the iron 6.4-keV 
line profile distortion, yields a/M = 0.989;^ |^^^ [79] (see 
also Ref [801). As a result, the kinetic energy of rotation 
fkin = Jr^ /2 storcd in the central engine turns out to 
be quite sufficient to explain the energy source of activity 
of compact objects. 

There are no particular problems with a regular mag- 
netic field, either. In young stars, the proper magnetic 
field Bq is measured directly and can be as high as 
lO'^ G [71]. At present, there are no direct observations 
of magnetic fields in radio pulsars, but they can be mea- 
sured in X-ray (accreting) pulsars, which are also neu- 
tron stars [11]. Therefore, nobody now doubts that the 
magnetic field of a neutron star can reach 10^^ G, and 
even extend up to 10^^ G in magnetars [81]. The situ- 
ation is worse with the magnetic fields of supermassive 
black holes. As is well known, a black hole cannot have 
a proper magnetic field, but the field can be generated in 
the surrounding accretion discs [82]. Unfortunately, so 
far there is no self-consistent theory of such generation, 
so we have to apply the estimate Bq ^ B^dd, where 

^-'^-^"^lo^) 

Let us keep in mind that such an estimate comes from 
the simple assumption that the energy density of the 
magnetic field is comparable to the total energy density 
in the accreting plasma yielding the Eddington luminos- 
ity (2). Clearly, estimate (7) represents rather an upper 
limit of the magnetic field near the black hole. In par- 
ticular, it does not take into account the contribution 
from the thermal pressure, which can be significant in 
gamma-ray burst sources. 

Finally, the problem of the 'electric wiring' can also 
be easily solved at first glance. Due to the presence of a 
strong magnetic field, in all cases the Larmor radius of 
particles tl = mcv / eB is always much smaller than the 
size R of the central engine. Therefore, one can consider 
with good accuracy that the electric current flows along 
the direction of the regular magnetic field. However, here 
we meet the problem of current closing, since particles 




separatrix 



braking neutron star 
torque K 



Figure 7: The structure of electric currents (contour ar- 
rows) near the polar caps of a neutron star. The Am- 
pere force related to the surface current Jg, produces 
the torque K braking the neutron star rotation. Above 
the acceleration region, the energy flux is predominantly 
transported by the Poynting vector (hatched arrows) . 



in the region of the load must move across the magnetic 
field. We shall necessarily discuss this point below. 

As an example. Fig. 7 illustrates how the braking 
occurs in an axisymmetric magnetosphere of radio pul- 
sars. Clearly, the total current fiowing out of the pul- 
sar surface must vanish; thus, there must necessarily be 
a reverse current in the magnetosphere to compensate 
for the loss of charges from the neutron star. As a re- 
sult, currents Jg closing the longitudinal currents in the 
magnetosphere must fiow over the pulsar surface. The 
ponderomotive action of these currents must brake the 
rotation of the radio pulsars [3, 56]. 

Thus, the problem of the magnitude of potential dif- 
ference is solved quite easily. But the problem of the 
load that determines the current / and, hence, the en- 
ergy losses, proved much more difficult. A long way had 
to be covered in order to solve it, and this, essentially, 
will be discussed in this review. Nevertheless, we shall go 
somewhat ahead and give here the preliminary estimates 
confirming the applicability of the discussed mechanism. 
As shown below, a good estimate of the electric current 
density is given by the expression 

JGJ = PGJC, (8) 

where 

Po..--^ (9) 

is the electric charge density that is needed for the elec- 
tric field in the rotating reference frame to vanish. For- 
mula (9) can be easily derived from relation (4). It was 
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first applied to the neutron star magnetosphere in the 
pioneering paper by P Goldreich and W H Julian [3], 
so the charge density (9) is usually called the Goldreich 
density. 

Clearly, the total electric current circulating in the 
magnetosphere of the central engine can be conveniently 
written out in the form 



^tot = iolcj- 



(10) 



Here iq is the dimensionlcss current, and Iqj = T^H^cpcj, 
i.e., for the case pcj ~ const wc obtain 



/gj = 



(11) 



Finally, Rq is again the size of the working area on 

the central engine surface. For black holes we can set 
i?o ~ i? = Tg, and for neutron stars (radio pulsars) it 
must be on the order of the radius of a polar cap from 
which magnetic field lines can go beyond the light cylin- 
der Ri, = c/Cl. Indeed, inside the closed magnetosphere, 
by virtue of the remarkable Ferraro isorotation law, the 
plasma starts rotating with the star as a solid body, and, 
hence, this region cannot work as a unipolar inductor. 
The working area will include only the region of open 
field lines, inside which the plasma rotational velocity 
can be different from that of the star. As a result, for 
the dipole magnetic field we obtain 



Ro 



Rl^) 



1/2 



(12) 



For relativistic strongly magnetized wind it is natural 
to assume that 

io « 1, (13) 

which corresponds to a free plasma outflow with the ve- 
locity c. As we shall see, this estimate is indeed correct. 
Therefore, the total energy losses can be estimated as 



tot 



nRo 



b^rIc. 



(14) 



In consequence of this, as shown in Table 1, the unipolar 
inductor model allows us to explain both the total energy 
release Wiot and the time of activity of compact sources, 
td — ^kin/W^tot- As mentioned above, for radio pulsars 

estimate (14) with account for relation (12) coincides to 
within an order of magnitude with the magnetodipole 
losses (3). 

For nonrelativistic outflows, estimate (13) is incorrect, 
and, as a detailed analysis shows, io ^ 1 [9]. For a 
sufiiciently rapid rotation with Cl > Q^t, where 

a. = ^^^^V''~10-s- (15) 



Ro \ Bf 







we have 



c / rip 



^in \ ^cr 



-2/3 



(16) 



and the dimensionless current for slow rotation limit 
takes the form 

io « — . (17) 

Here pin is the density of the outflowing matter near the 
surface of the star, and t^in is the characteristic velocity 

of the outflow along the jet axis. As a result, the total 
energy losses for rapidly rotating stars can be expressed 
through the directly observed quantities: 



(18) 



i.e., through the total magnetic flux ^tot = t^RoBo, the 
rotational angular velocity fi, and the mass loss rate in 
the jet M. For the parameters typical in young stars we 
have 



10 



,36 



10" cm 



8/3 



P 

lOS 
M 



4/3 



Bu 



1/3 



4/3 



10-9 Mq yr- 



ergs 



(19) 



It is seen that this value is indeed close to the energy 
losses from young stellar objects. Thus, the unipolar 
inductor model allows us to explain the main jet charac- 
teristics for nonrelativistic sources, too. 

Interestingly, the knowledge of the total energy losses 
Wtot immediately allows the total longitudinal electric 
current the circulating in the magnetosphere to be esti- 
mated. Indeed, by comparing expressions (11) and (14), 
we straightforwardly obtain 



(20) 



The characteristic amplitudes of currents are also col- 
lated in Table 1. 

3.2 Grad-Shafranov equation method 

The Grad-Shafranov equation method lies at the heart 
of the analytical theory which, in our opinion, is able 
to quite successfully describe the main properties of ac- 
tive compact astrophysical sources. Simply speaking, 
this approach describes axisymmetric stationary flows 
in the framework of ideal magnetohydrodynamics. This 
approximation is based on the assumption of a high con- 
ductivity of the plasma that fills the magnetosphere of 
the central engine (the high energy release guarantees 
a high degree of ionization of matter, and the effective 
production of electron-positron pairs in the vicinity of 
black holes). Moreover, most of the sources discussed 
above (except for radio pulsars) can be considered to a 
good approximation as axisymmetric and stationary. 

The attractiveness of this approach is related to the 
fact that there arc quite a lot of integrals of motion in 
stationary ideal magnetohydrodynamics, i.e., quantities 
which are conserved along particle trajectories. This im- 
mediately provides us with important information with- 
out complicated calculations. Indeed, to determine the 
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Table 1: Parameters of the central engine: AGN — active galactic nucleus, GRB — gamma-ray burst, /iQSO — 
microquasar, PSR — radio pulsar, msPSR — millisecond radio pulsar, and YSO — young stellar object. 





AGN 


GRB 


HQSO 


PSR 


msPSR 


YSO 


Mass M in M© 


IC-IO^ 


- 10 


- 10 


w 1.4 


w 1.4 


- 10 


Radius R, cm 




- 10" 


- 10" 


- 10" 


- 10" 


^10^^ 


Working radius Rq 


- R 


- R 


- R 


{VlR/c)^/'^R 


[VLR/cyr^R 


- R 


Period P 


10-10^ s 


^ 1 ms 


^ 1 ms 


- 1 s 


1.39-10 ms 


1 10 d 


nR/c 


> 0.1 


> 0.1 


> 0.1 


- 10-* 


^10-^ 


- 10-" 


Magnetic field Bq, G 


lO'^-lO'' 


- 10^=' 


^ 10^*^ 


^10^^ 


^ 10« 


- 10^ 


Energy storage fkin, erg 


lO^^-lO"^ 


- 10^^ 


- 10^^ 


10**~10*" 


^10^^ 


- lO'^* 


Dimensionless current iq 


1 


1 


1 


1 


1 


^ Cj Win 


Power Wtot, erg s"^ 


10^^-10*=' 




^ 10^« 






- 10^" 


Lifetime td, yr 


^10'^ 


- 10-" 


- lO'^ 


10"-10^ 


IG^-IG'J 


- 10* 


Current /, CGSE 




-10^^ 


- 10^" 


10^^-10^^ 


- 10^^ 


- 10^" 



T(r,0)=0 




Figure 8: 
const. 



Axisymmetric magnetic surfaces ^'(r. 



height a throwing stone reaches it is not necessary to 
solve equations of its motion: it is sufficient to apply the 
energy conservation law. 

In the axisymmetric case, as illustrated in Fig. 8, the 
magnetic field vectors must lie on the magnetic surfaces 
which can be easily parametrized using the magnetic flux 
function ^(r., 9) that determines the magnetic field 



^ V^- X 

£> = 

27rn7 



21 



(21) 



Here zu = r sm6 is the distance from the rotational axis, 
and the numerical coefficient in the first term is chosen 
such that the function '^{r,6) indeed coincides with the 
magnetic flux passing through a circle r, 9, < if < 27r. 
As for the quantity I{r,9), it represents the total elec- 
tric current flowing through the same circle. It is easy 
to check that the following important properties are sat- 
isfled. 



1. At all times d^P = B • dS (dS is the surface ele- 
ment). Therefore, the function 'ii{r,9) indeed bears 
the sense of the magnetic flux. 

2. Since the poloidal part of the magnetic field in for- 
mula (21) can be written out as (27r)-^V^'x Viy9, the 
condition V-B = is automatically satisfied. Thus, 
three components of the magnetic field are com- 
pletely determined by two scalar functions ^'(r, 9) 
and /(r, 9). 

3. For the same reason it is clear that the condition 
B • V^* = will be satisfied for axisymmetric flows. 
Therefore, the lines 5'(r, 6*) — const deflne the form 
of magnetic surfaces. As a result, the integrals of 
motion should depend on only one scalar function, 
*(r,6'). 

Let us now understand which integrals of motion ap- 
pear in the case of axisymmetric stationary flows. Inci- 
dentally, the very structure of the Grad-Shafranov equa- 
tion method will become clear. For simplicity, let us 
consider first a purely hydrodynamic flow. In this case, 
in analogy with relation (21), it is necessary to introduce 
the function $(r, 9) of hydrodynamic flux deflned as 



X e, 



27rr sin 9 



(22) 



where hereinafter the subscript 'p' will correspond to 
the poloidal [i.e., lying in the (r, 9) plane] components of 
vectors. 

In hydrodynamics, there are flve scalar equations (the 
mass continuity equation, three components of the mo- 
mentum Euler equation, and the energy conservation 
equation) for five unknown quantities — three velocity 
components, and two thermodynamic functions. How- 
ever, due to the axial symmetry, stationarity, and ideal- 
ity of the flow three of five equations can be represented 
in the form (vV)!*^'-' — 0, which means that integrals I*^*^ 
must be constant on the surfaces ^{r, 9) — const. As is 
well known, these integrals include the energy (Bernoulli 
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integral) E^, the specific angular momentum in, and the 
entropy s: 



in 



X„($) = w^rsin^, 
s($). 



(23) 

(24) 
(25) 



Here, <^g is the gravitational potential, w{p, s) is the 
specific enthalpy, and the subscript 'n' corresponds to 
nonrelativistic quantities, while Bernoulli integral -E'n('^) 
corresponds to the projection of the Eulcr equation (i.e., 
the equation of the force balance) onto the direction 
along the poloidal velocity Vp, and the angular momen- 
tum ivn($) represents the projection onto the unit vector 
e^. The remaining two first-order equations can be re- 
duced to one second-order equation for the flux function 
<i>(r, (?). It is clear that this equation will describe the 
force balance in the direction perpendicular to surfaces 
$(r, 0) = const. In the compact form, it can be written 
out as 



— zu^Wk 



d$ 



T ds 



d<I> mp d<I> 



0. 



(26) 



It should be noted that we deliberately defined the 
enthalpy w as a function of density p and entropy s. The 
point is that Bernoulli equation with the use of definition 
(22) can be recast into the form 

One can see that Bernoulli equation written in this form, 
in addition to the integrals of motion and the flux frmc- 
tion <^{r,9), contains only the density p. Consequently, 
it indirectly deflnes the density p though the fiux func- 
tion 4* and integrals of motion: 



p(V<I>;£;n,L„,s;r,( 



(28) 



This implies that after substituting Eqn (28) into Eqn 
(26), the latter will contain only one unknown fiux func- 
tion $(r, ^) and three integrals of motion depending on 
it. 

Such an equation can also be written out in the frame- 
work of ideal magnetohydrodynamics, too. In this case, 
however, not three but five integrals of motion exist. 
Two additional integrals come from the freezing-in con- 
dition 

E-|-vxB/c = 0. (29) 

Indeed, from condition (29) follows that the electric 
field is perpendicular to the magnetic field. In the 
axisymmetric case this means that magnetic surfaces 
^{r,9) = const will be equipotential ones. This condi- 
tion can be conveniently rewritten in the form 



Of = Of(*), 



(30) 



where the scalar quantity Of determines the electric field 
according to the definition 

Of 



E 



27rc 



(31) 



This is related to the fact that: 

• in the axisymmetric case {d/dt = 0) Maxwell equa- 
tion V X E = leads to the condition = 0, 

• the freczing-in condition yields = 0, 

• the definition [31] together with Maxwell equation 
V X E = leads to the condition VOp x V'J' = 0, 

where relation (30) comes from. 

Function Of introduced in this way bears the meaning of 
the angular velocity of particles (more precisely, the mo- 
tion of particles is the sum of rotation with the angular 
velocity Of and sliding along the magnetic field). Con- 
dition (30) represents the Ferraro isorotation law [83], 
according to which the angular velocity of particle mo- 
tion relative to the magnetic field must be constant on 
axisymmetric magnetic surfaces. 

On the other hand, the freezing-in plasma condition 
implies that plasma velocity vectors v have also to lie on 
the magnetic surfaces, i.e., the flux of matter does not 
intersect the boundaries of the magnetic surfaces. This 
means that the particle flux function ^{r,d) must be a 
function of the magnetic fiux \['(r, 9). This fact allows us 
to introduce one more integral of motion 



d<I> 



(32) 



which, as evidenced by the foregoing, bears the sense of 
the ratio of the particle flux to the magnetic field flux. 
Correspondingly, the poloidal velocity of matter can be 
written as 



(33) 



As for the energy and angular momentum integrals 
(which in the nonrelativistic case are usually considered 
as functions of the particle flux), they now take the form 



Enm = 

in($) = 



Of/ 



+ 



+ W + (fg, 



2wrinC 
h Vu,rsm6 

2'KrjnC 



(34) 
(35) 



respectively. The entropy s(\E') is ones again the one 
more (fifth) invariant. It is clear that both particles and 
electromagnetic field contribute to the energy and angu- 
lar momentum, and, as can be easily checked, the term 
OF//27r?7nC corresponds simply to the Poynting vector 
flux. 

Next, Bernoulli equation (34) can now be rewritten in 
the form 



(1-7W2)2 



2vD^^Y 



(36) 
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where 



P 



(37) 



The quantity is the square of the Mach number of 
the poloidal velocity v-p relative to the poloidal compo- 
nent of the Alfven velocity: 



Br. 



(38) 



i.e., M.^ — It should be recalled that the specific 

enthalpy w in equation (36) must be considered as a 
function of entropy s, as well as of the Mach number 
Ai^ and the integral r^j,. The corresponding relationship 
has the form 



X2 



dP_ 

ds 



T 



Vs. 

(39) 

In consequence, as in the hydrodynamic limit, Bernoulli 
equation allows one to determine, albeit indirectly, the 
quantity via the magnetic flux ^(r. 0) and five inte- 
grals of motion: 

= M^Vi';E^,L^,s,7]^,QF;r,e). (40) 

As for the projection of the force balance onto the di- 
rection perpendicular to the magnetic surfaces, it can be 
written in the forrrQ 



1 



1 - 

2{E^ - w 

. J_d77n J 

ryn d^P n 



df^F 
"d^ 



d^n 
d* 



1 M^L^-^^uP- AL^ 



^2 \~M? 



d* 



^2(1- 7^2) 



:S-»' 



(41) 



Now the structure of the treatment considered here 
becomes clear. Equation (41) jointly with Bernoulli 
equation (36) determines the value of the magnetic flux 
^'(r, 0). Then, again using Bernoulli equation, one can 
determine the value of the Mach number M. at each 
point. After that it turned out that the number of inte- 
grals of motion is enough for determining all other quan- 
tities from simple algebraic equations. For example, one 
arrives at [841 



/ 



1-7W2 ' 

1 »Fn72 _ l^^j^l 
w 1 - A42 ' 



(42) 
(43) 



and, respectively, p = Anrj^/M.'^ . This is essentially the 
main attractiveness of the approach we discuss. Some- 
times, as we shall see, the key properties can be directly 



^ Unfortunately, in monograph | 
were discarded. 



-</3g in Eqn (4.102) 



obtained from algebraic relations. In other words, by 
making sufficiently reasonable assumptions about the 
structure of the fiow [i.e., about flux function '^{r.O)], 
it is possible not to solve equation (41) at all and to 
analyze only algebraic, albeit indirect, relations. 

The full version of the nonrelativistic equation con- 
taining all flve invariants was flrst formulated by L S 
Solov'ev in 1963 in the third volume of the Reviews of 
Plasma Physics [85]. Being virtually unknown for as- 
trophysicists, this equation was later reformulated anew 
several times [86-88]. For this reason, in particular, to 
date there has been no unique system of notations, so 
sometimes it is difficult to compare the results of dif- 
ferent studies. In the literature, equations of this type 
are commonly called Grad-Shafranov equations, which 
were formulated at the end of the 1950s in relation to 
controlled thermonuclear fusion [89, 90], although the 
hydrodynamic version of this equation was known even 
earlier (see, for example, Ref. [91]). Similar equations 
going back to the classical Tricomi equation, were dis- 
cussed as early as the beginning of the twentieth century 
in the context of transonic hydrodynamic flows [92, 93]. 

For simplicity, we have written out above equations 
only for the nonrelativistic case. However, it was not too 
difficult to obtain the corresponding equations both for 
the relativistic case [94] and for flows in the vicinities of 
nonrotating [95] and rotating [96, 97] black holes, since 
the Kerr metric is axisymmetric and stationary. It is 
these relativistic equations that we shall discuss below. 
In the main text we shall try to formulate sufficiently 
simple asymptotic expressions by focusing on the qual- 
itative description of the flow properties. Sufficiently 
lengthy full equations are presented in the Appendix. 
Here, we shall restrict ourselves by writing out addition- 
ally the integrals of motion for relativistic flows in flat 
space. 

Clearly, magnetic surfaces remain equipotential in the 
relativistic case, too. Thus, the angular frequency fip 
in deflnition (31) remains the integral of motion. As for 
the integrals of energy E and the z-component of the 
angular momentum L, now they should be written out 
afl 

E = E{-^) = ^+-fp,i^c^, (44) 

L = i(^') = + p,j]wuu,c. (45) 

Here u is the spatial part of the four-velocity vector 
(7 = Vu2 + 1 is the Lorentz factor) , and 



p 



nipC + nipW + 



(46) 



is the relativistic enthalpy including the rest mass of par- 
ticles. Finally, the relativistic integral of motion 77 is now 
determined from the condition 



P' 



(47) 



^To avoid misunderstanding, from now on the electric field E 
will always be boldfaced. 
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where n is the particle number density, and hereinafter 
all thermodynamic functions will be defined in the co- 
moving reference frame. Thus, for relativistic flows 
where |up| « 7 we simply have 



rj - 



^(lab) 

'b7' 



(48) 



It should be noted that the relativistic and nonrelativis- 
tic integrals of motion have different dimensions, since 
the relativistic integrals are normalized not on the unit 
matter flux d<i> but on the unit magnetic flux dvP. 

It is evident that again both the energy flux and the 
angular momentum flux include the contributions from 
the electromagnetic field and particles, with the electro- 
magnetic contribution [accurate up to an additional fac- 
tor 77(5')] fully coinciding with that obtained in the non- 
relativistic limit. In the general magnetohydrodynamic 
case, the total energy and angular momentum losses, 
Wtot and -ftTtot, will be determined by the relationships 



1 



£;(*)d«'. 



L(*)d*, 



(49) 
(50) 



respectively. 

3.3 Supersonic flows 

3.3.1. The model. In order to clearly understand the 
main features of the model considered here, it is conve- 
nient from the very beginning to analyze a sufficiently 
simple geometry of magnetic surfaces, and to formu- 
late basic parameters characterizing the flow. Figure 
9 demonstrates the simplestsplit monopole model of the 
magnetized wind [33], which during many years served 
as the 'hydrogen atom' for all researchers who studied 
the nature of the activity of galactic nuclei, gamma-ray 
bursts, and microquasars. Notably, most analytical re- 
sults were obtained exactly for such flows. 

It is assumed in the framework of this model that 
the 'central engine' involves a compact object (neutron 
star or black hole) and an accretion disc which separates 
the converging and diverging magnetic field fluxes. The 
accretion disc here is nedeed both to separate the op- 
positely directed magnetic field fluxes and, in the case 
of a black-hole magnetosphere, to generate the regu- 
lar poloidal magnetic fleld (it will be produced by the 
toroidal currents flowing in the disc). In the absence of 
the accretion disc, a black hole, as is well known, cannot 
have the proper magnetic field (the so-called 'no-hair 
theorem' [18]). In addition, poloidal currents will also 
fiow in the disc, closing the bulk currents flowing in the 
magnetosphere. Notice that a similar conflguration with 
a disc separating magnetic field fluxes far from the neu- 
tron star also emerges in many models of radio pulsars 
[98 104], because it is natural to assume that at large 
distances the flow becomes quasispherical. 




Figure 9: The structure of electromagnetic flelds for 
the split monopole magnetic fleld near a slowly rotating 
black hole [33] . Currents flowing in the highly conduct- 
ing disc in the equatorial plane provide both the poloidal 
magnetic fleld jump and the closure of bulk currents flow- 
ing out the upper and bottom hemispheres. 



Next, it is important that the crossed flelds Eg and B,^ 
form the electromagnetic energy flux (the Foynting vec- 
tor flux), which is directed along the magnetic surfaces. 
It should be stressed that this energy is transferred at 
the zero frequency, so the electromagnetic field that car- 
ries the energy is not an electromagnetic wave in the 
usual sense. The electromagnetic energy flux therefore 
appears only due to the longitudinal current generating 
the toroidal magnetic field; in the absence of particles, 
such energy release becomes impossible. The plasma also 
moves along the magnetic surfaces, so the sum of their 
energy fluxes is the integral of motion. On the other 
hand, the longitudinal current / is not the integral of 
motion, so the MHD approximation we are considering 
allows, in principle, describing the current closure phe- 
nomenon. However, the magnetic surfaces here remain 
equipotential. Therefore, such flows can carry high elec- 
tric voltages over large distances from the central engine. 
This should always be borne in mind when discussing the 
interaction of a magnetized wind with the surrounding 
medium. For example, this effect must be taken into 
account in close binary systems with radio pulsars. 

As for the main dimensionless quantities characteriz- 
ing the flow, they include the magnetization parameter 
(T, the particle production multiplicity A, and the com- 
pactness parameter la- The magnetization parameter a 
shows by how much the electromagnetic energy flux near 
the central engine can exceed the particle energy flux. 
Thus, as seen from the definition of the energy integral 
(44), it can be more conveniently defined for relativistic 
flows as 

^ ^ (51) 



where the maximal value is chosen for all magnetic sur- 
faces. As a result, the value of a corresponds to the max- 
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imal Lorentz factor of the plasma that can be reached 
in the case where all the electromagnetic field energy is 
transferred to the particles. In other words, a is the max- 
imum Lorentz factor that can be achieved in the mag- 
netized wind. Of course, here the mean hydrodynamic 
energy of the plasma fiowing out is assumed. In particu- 
lar, for the split monopole magnetic field (for which this 
quantity was first introduced by F C Michel [4] in 1969), 
we obtain 

02,Tf 

(52) 



a = 



Correspondingly, for the nonrelativistic flow it is conve 
nient to use the quantity 



(53) 



where f in is tlic velocity of matter flow along the jet axis. 
It is easy to check that the strong magnetization condi- 
tion r2 > Qcr (15) coincides with the condition an > 1. 

Recall that we are mainly interested in strongly mag- 
netized flows, i.e., flows in which the main energy flux 
near the central engine is due to the Poynting vector 
flux fl-pl /2tt.. In the opposite case, the flow will be only 
slightly different from the hydrodynamic outflow. Using 
definition (51), this condition can be rewritten as 



7iii < f^, 



(54) 



where 7in is the injection Lorentz factor. As we shall 
see, the magnetization parameter a is the key parameter 
determining the basic features of the flow. 

Next, to find the ejected plasma density, it is conve- 
nient to introduce the dimensionless particle production 
multiplicity A 

(lab) 

A , (55) 

whereriGJ = |PGj|/e. Such a deflnition is connected with 
the fact that, as we shall show below, both in pulsar mag- 
netospheres and in black hole magnetospheres the dens- 
est is the secondary electron-positron plasma generated 
either due to conversion of hard gamma-quanta in the 
magnetic fleld or due to gamma-ray collisions with ther- 
mal photons [33, 105]. However, hard gamma-quanta 
must be emitted in both cases by primary particles whose 
density is supposed to be close to the Goldreich density. 
The convenience of the particle production multiplicity 
A also stems from the fact that the magnetization pa- 
rameter a can be rewritten with its help in the form 



<T = 



4AmeC^ 



1 /Wu 



Wa 



1/2 



(56) 



where Wa = m^e' /e^ ^ IQ^'' erg s~^. Thus, the knowl- 
edge of two of three quantities Wtot, and A, allows the 
determination of the third one. 

Finally, the compactness parameter 



is in fact the optical depth for Thomson cross section 
fTT at a distance R from a source with the total lumi- 
nosity Ltot- Below, it will be important for us that the 
parameter provide an upper limit of particle energy in 
the acceleration region. On the other hand, a large la is 
necessary for effective particle production. 

3.3.2. Singular surfaces. Singular surfaces rep- 
resent the most important structural element of flows. 
As we shall see, it is the analysis of the conditions of 
the smooth passing of a flow through singular surfaces 
that allows sometimes rather general relationships to be 
obtained without solving the Grad-Shafranov equation 
itself. Notice from the very beginning that, for simplic- 
ity, below we shall only analyze the case of cold flows. 
The point here is that, at large distances thermal effects 
are insigniflcant for the polytropic index F > 1 (pressure 
P (X nF). This conclusion can readily be obtained from 
both the Grad-Shafranov equation itself and Bernoulli 
equation. Indeed, from the analysis, for instance, of non- 
relativistic equations (36) and (41) it follows that both 
the enthalpy w = Cg/(F — 1) oc n^~^ in Eqn (36) and 
the temperature T oc nF~^ in Eqn (41) decrease with 
the distance from the compact source, since for any di- 
vergent outflow the particle number density n — )■ for 
r — > oo. Therefore, the contribution from a final temper- 
ature (enthalpy, entropy) compared to the total energy 
E and its derivative di?/d^' can be neglected at large 
distances. Thus, it becomes clear why in the analysis 
of relativistic fiows the final temperature effects (and, in 
particular, critical conditions on the slow magnetosonic 
surface) are usually neglected. On the other hand, the 
pressure can be significant for cylindrical flows, i.e., for 
flows in which the density does not decrease with dis- 
tance from compact object [9,106]. 

The first natural scale that emerges in the theory of 
relativistic winds is the light cylinder 



Ri, = 



(58) 



(57) 



i.e., the axial distance at which solid-body rotation to- 
gether with the central object becomes impossible. It is 
easy to show that the light cylinder is the scale where: 

1. the magnitude of the electric field becomes compa- 
rable to that of the poloidal magnetic field; 

2. toroidal electric currents flowing in the magneto- 
sphere start perturbing the poloidal manetic field of 
the central engine; 

3. the magnitude of the toroidal magnetic field pro- 
duced by the longitudinal Goldreich current be- 
comes comparable with that of the poloidal mag- 
netic field. 

It follows from the first statement above and definitions 
(21) and (31) that beyond the light cylinder the elec- 
tric field becomes stronger than the poloidal magnetic 
field. In particular, the poloidal magnetic field will de- 
crease as and the electric field as r~^ for the spit 
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monopole outflow shown in Fig. 9. On the other hand, 
freezing-in condition (29) requires that the magnetic field 
be stronger than the electric field. This can be possible 
only when a strong enough longitudinal electric current 
is flowing in the magnetosphere, because the toroidal 
magnetic field also decreases as according to Eqn 
(21). 

We can therefore conclude that the question as to 
whether or not a smooth relativistic MHD outfiow 
(|E| < |B|) exists beyond the hght cylinder is also di- 
rectly related to the question of the magnitude of the 
longitudinal current circulating in the magnetosphere of 
a compact object. Then, for currents below some critical 
value, a so-called light surface is bound to appear in the 
magnetosphere, on which the electric field matches the 
magnetic field (|E| — |B|), and hence the approximation 
considered here itself becomes invalid. Calculations [107, 
108] showed that the closure of currents occurs near this 
surface in the region with the thickness Sr ^ Ri^/\, and 
particles are effectively accelerated there up to energies 
of 7 ^ cr. 

If the longitudinal currents are sufficiently high, the 
smooth MHD outflow can exist beyond the light cylinder 
as well. The electric field there will be almost equal 
to the magnetic field. Indeed, as directly follows from 
the relativistic Bernoulli equation (^.12), in the limit 
w ^ i?L we obtain simplj(3 



lEl' 



^2- 



Since, as can be easily checked, Bp « B^/xr 
Xr = Q,vj / c, we can always apply the estimate 



Bt 



f>2 

|Ep< ^. 

7 



(59) 



where 



(60) 



As a result, the radial drift motion in the crossed elec- 
tromagnetic fields dominates in a strongly magnetized 
relativistic outflow beyond the light cylinder. Indeed, as 
can be easily verified, the Lorentz factor entering into 
Eqn (59) satisfies the condition 7^^ = 1 — U^j., where 

E X B 



Udr = C- 



52 



(61) 



In other words, the velocity parallel to the magnetic field 
does not contribute at all to the value of the Lorentz 
factor [109]. 

The fast magnetosonic surface is another important 
surface of the magnetized flows. It is fully equivalent to 
the sonic surface in the ideal hydrodynamics. Indeed, 
Bernoulli equation (27) is well known to have a singu- 
larity on the sonic surface. For example, the logarithmic 
derivative of density determined from Eqn (27) is written 
for a spherically symmetric flow in the form 



2v'- 



GM 



r dp 

^'1 " ~T' = 2 

p ar Cg 



2 - 



GM 



rv^ 



V 



4, (e. 



^This expression corrects Eqn (4.144) from monograph [9]. 



It is obvious that derivative (62) has a singularity 
when the velocity of matter equals the speed of sound: 
V = Cs = Ci. {D = 0). This means that in order to cross 
the sonic surface r ^ r^ smoothly, the additional condi- 
tion 

, , GM , , 

iVn ==2 ^=0 63 

r^ci 

must be satisfied. As a result, the additional critical con- 
dition (63) fixes the accretion (ejection) rate of matter 

[14]- 

The fast magnetosonic surface plays a similar role. 
But now it determines not the accretion or ejection rate, 
but the magnitude of the longitudinal current / (more 
precisely, the integral V). It is this critical condition on 
this surface that shows us that in the relativistic case the 
longitudinal current / near the central engine must be 
close to the Goldreich current /qj. On the other hand, 
as we have already noted, for a nonrelativistic outflow 
io ^ 1, and the conditions of the smooth crossing of 
singular surfaces lead to relations (15)-(17) used above. 
Their derivation, however, is rather cumbersome, and 
we shall omit it here. It should only be emphasized 
that they can be obtained directly from an analysis of 
Bernoulli equation. The point is that the sonic surface is 
the X-point on the (distance r-velocity v) plane. That 
is, it is the point of crossing the roots of the algebraic 
Bernoulli equation. The condition of coincidence of roots 
of the algebraic equation puts certain bounds on the co- 
efficients of the equation itself, which enables the magni- 
tude of the longitudinal current to be estimated. Notice 
that expressions for the current formulated above were 
exactly obtained in Refs [110,111] for the simplest split 
monopole geometry shown in Fig. 9. 

In a similar way, the following theorem can be proved: 
In a relativistic outflow near the outer fast magnetosonic 
surface, the energy of the particles reaches the values 



( ^ 



7 



\pr]c? 

7im 



1/3 



o^l\ 7in«^^/^ (64) 



7in > O- 



1/3 



(65) 



Thus, the fraction of energy carried by particles in the 
vicinity of the fast magnetosonic surface is a small frac- 
tion (~ (T~2/^^ of the electromagnetic energy flux for 
strongly magnetized outflows a 3> 7in- The surface it- 
self is located at the distance of 



( E 



\ firjc^ 
E 



1/3 



Rl - a^/^RL, 7in « (7^/^(66) 



1/3 



1/2 



Rl 



1/2 



Rl, 



(67) 



7in > CT 



1/3 



(the first relation holds true not too close to the rota- 
tional axis). Interestingly, expression (66) is valid for 
both relativistic and nonrelativistic flows since it does 
not include in fact the speed of light c. 
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Finally, the singularity a.t A = 1 — = that ap- 
peared in nonrelativistic equations (42), (43) suggests 
that the Alfvenic surface must also play an important 
role in the structure of magnetized flows. Its location 
for nonrelativistic flows can easily be estimated from the 
numerator of relation (42): 



-^2 



In this case, the Alfvenic surface turns out to be located 
close to the fast magnetosonic surface. 

As to a relativistic flow (and a flat space), the corre- 
sponding condition should be written differently: 



A^l 



(69) 



On the other hand, it is easy to check that the parame- 
ter q = Wpart/Wem (i-G., the particlc-to-electromagnetic 
energy flux ratio) can be presented in the form 



q = 



(70) 



Thus, in the region where the energy flux is Poynting- 
dominated {q <C 1), the following condition should be 
satisfied: 



(71) 



Consequently, the Alfvenic surface for such flows is lo- 
cated near the light cylinder. Hence (except for the 
polar region where both the Alfvenic and fast magne- 
tosonic surfaces are close to each other), the fast magne- 
tosonic surface for strongly magnetized flows is located 
cr^/"^ times further from the central engine as compared 
with the Alfvenic surface. 

Let us remember that the Alfvenic surface in the rel- 
ativistic case determines the scale on which the toroidal 
magnetic field becomes comparable in magnitude to the 
poloidal field. It is easy to check that for rapid rota- 
tion 57 > ilcr (15) a similar situation holds for nonrel- 
ativistic flows, too. As regards the electric field, it is 
always weaker than the magnetic field in the nonrela- 
tivistic case. Notably, that is why the light surface can- 
not appear in the nonrelativistic case. 

The Alfvenic surface represents a higher-order singu- 
larity as against the fast magnetosonic surface. There- 
fore, relations (42) and (43) do not put any constraint 
on the integrals of motion and only determine the lo- 
cation of the Alfvenic surface and the magnetic field 
structure. In this event, however, particles can intersect 
the Alfvenic surface only in one direction. For example, 
when the central source loses its rotational energy, cross- 
ing the Alfvenic surface is possible only in the direction 
outward from the compact object. When the energy flux 
is directed toward the central engine (for instance, if it 
is spun up by the accreting material) , the flow must also 
be directed toward the central engine. Certainly, this 




Figure 10: The location of the Alfvenic (A) and fast 
magnetosonic (F) surfaces near the black hole horizon. 
The dashed line shows the Alfvenic surface in the force- 
free approximation (i.e., the 'light cylinder'), while the 
dotted line indicates the ergosphere surface. Here, a is 
the gravitational redshift (see the Appendix). 



statement is invalid in the region of an accretion disc 
where viscosity cannot be neglected. 

This statement can easily be proved in the relativistic 
case by recalling that the motion of particles is the sum 
of the drift motion in the crossed fields and the motion 
along the magnetic field. The condition v < c that lim- 
its the longitudinal velocity puts bounds on the radial 
velocity of matter. This is related to the fact that the 
drift velocity itself becomes close to the speed of light on 
the Alfvenic surface. However, accretion of matter with 
positive energy release cannot be realized in the nonrel- 
ativistic case, either. In this event, the interaction of the 
supersonic accretion flow with the rotating magnetized 
central body would take place. A shock wave is known 
to be formed in such an interaction [112]. 

In conclusion, we should comment on the features of 
a black hole magnetosphere. As seen from exact ex- 
pressions for the Alfvenic A (A. 13) and sonic D (^.20) 
factors presented in the Appendix, a second family of 
singular surfaces inevitably emerges near the black hole 
horizon. Here, as shown in Fig. 10, the infalling mat- 
ter, as in the case of the outflow, must first cross the 
Alfvenic surface and only then intersect the fast magne- 
tosonic surface (recall that thermal effects are not dis- 
cussed here). This is related to the fact that the strong 
gravitational field of the black hole forces the matter to 
approach the event horizon. 

The appearance of the second family of singular sur- 
faces leads to new important properties. First of all, the 
matter can intersect the inner Alfvenic surface only in 
the direction towards the black hole horizon. But this 
means that if the central engine loses rotational energy 
(and hence particles can cross the outer Alfvenic surface 
only in the direction away from the compact object), 
the plasma is bound to be generated in the magnetic 
field lines 'anchored' to the black hole horizon. Only in 
this case can electric currents appear in the black hole 
magnetosphere, which are necessary, as we have seen, 
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to explain the observed energy release. In turn, the ap- 
pearance of one more critical condition on the inner fast 
magnetosonic surface proved to be sufficient to deter- 
mine the angular velocity ftp. Here, ilp indeed must be 
close to riH/2 [r^H = '^(''g) is the angular velocity of the 
black hole; see the Appendix], as was understood as early 
as the Blandford and Znajek paper [33]. This problem 
has been possible to solve exactly for slow rotation in 
the split monopole magnetic field [113]. 

3.3.3. The problem setting. Before considering 
the main results which were obtained using the analyt- 
ical theory, it is necessary to discuss the formulation of 
the problem. The point is that we shall primarily be 
interested in transonic flows, i.e., those which are sub- 
sonic near the compact object and supersonic in the wind 
region. Indeed, as we shall see, the distance from a cen- 
tral engine to singular surfaces in all compact objects 
are much smaller than even the transverse size of the 
collimated outflows. The difficulty here is that the di- 
rect problem setup itself in the framework of the Grad- 
Shafranov equation method turns out to be nontrivial. 
For example, the second-order equation describing the 
two-dimensional flow in the hydrodynamic limit, when 
only three integrals of motion are available, requires 
four boundary conditions to be imposed in the transonic 
regime. The flfth condition is the critical condition on 
the sonic surface. This means that on some surface, for 
example, two thermodynamic functions and two velocity 
components must be specified. We stress that here only 
flows depending on two variables are considered. As dis- 
cussed in detail in monograph [9] , the well-known spher- 
ically symmetric flows (the Bondi accretion, the Parker 
ejection) are degenerate, since the structure of the flow 
itself is specified in them. In the general case of spheri- 
cal accretion, a nonstationary solution with a shock wave 
appears (see, for example, Ref. [114]). 

However, to determine Bernoulli integral, which is nat- 
urally needed in solving the equilibrium equation, we 
should specify all three components of the velocity, which 
is impossible since the third velocity component itself 
should be found from the solution. In the general case, 
one should set 

b = 2 + i-s' (72) 

boundary conditions, where i is the number of invari- 
ants, and s' is the number of singular surfaces (and for 
the magnetic field lines threading the black hole hori- 
zon this number is doubled, i.e., separately for ejecting 
and accreting plasma). Such an internal inconsistency 
of this approach in the general case does not allow us to 
solve direct problems, namely, to determine the structure 
of the flow in some region using the given physical pa- 
rameters on its boundary. Therefore, it is not astonish- 
ing that most researchers primarily interested in astro- 
physical applications already in the middle of the 1990s 
started addressing a totally different class of equations, 
namely those covering time relaxation problems, which 
can only be solved numerically [79, 115-119]. However, 
only in the last several years has signiflcant progress here 
been achieved [103, 109, 120-127], which, among other 



things, has conflrmed many analytical results obtained 
before. 

It should be noted that this problem does not ap- 
pear in both subsonic and supersonic cases. For these 
flows, all necessary integrals of motion must be deter- 
mined from the boundary conditions. In particular, the 
boundary conditions will determine in the subsonic case 
the longitudinal current /, too. This exactly corresponds 
to the unipolar inductor model, where the current (and 
hence the energy losses) is determined by the external 
load. Unfortunately, this ideology has also spread into 
the theory of magnetized winds. This is related to the 
fact that many results were obtained in the 1970s- 1980s 
using the force- free approximation [i.e., when cr — > oo 
and masses of particles can be neglected] . In this approx- 
imation, the Grad-Shafranov equation becomes elliptical 
and, hence, the flow structure must depend on conditions 
at the external boundary. But the theory of pulsar mag- 
netosphcres so far has been constructed in the force-free 
approximation. The class of subsonic flows also includes 
the so-called 'magnetic tower' [128, 129]. As this ques- 
tion is highly important, we shall discuss it below in more 
detail. 

Correspondingly, the Blandford-Znajek model was 
also constructed in the force-free approximation, which 
required the boundary condition to be set on the black 
hole horizon. Since the electromagnetic wave (like other 
material bodies) can propagate near the horizon only 
normal to the horizon toward the black hole, the bound- 
ary condition in fact is equivalent to the Leontovich 
boundary condition in radio physics [130]. However, it 
is usually obtained by requiring the finiteness of fields 
in the freely falling observer's frame of reference, which 
yields B^{rg) = —Eg{rg). This condition, as is well 
known, can be rewritten in the form of the Ohm law 
for the formally introduced 'surface current' [131] 

J(rg) = J^E(rg), (73) 

which corresponds to the universal 'internal' resistance 
of the battery, TZ = A-k/c = 377 O. In another form, this 
boundary condition can be rewritten as 



47r/(*) = [Oh - Of(*)] sin 6I 




Here, we have utilized definitions {A.7) and (A. 8) and, 
for simplicity, written the equality B^p{rg) = —Effifg) 
for slow rotation. Thus, it is not surprising that in the 
framework of this approximation the mechanism of en- 
ergy loss by a black hole was connected, in analogy with 
the unipolar inductor, with the Ampere force acting on 
the black hole horizon from the side of the surface cur- 
rent [131]. 

Only much later was it understood that the force-free 
approximation gives inaccurate and sometimes erroneous 
results. This is related to the fact that in the force-free 
approximation, i.e., when particles are assumed to be 
massless, the flow always remains subsonic. Under this 
assumption, the fast magnetosonic surface on which the 
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Figure 11: The structure of a magnetic field behind a 
switching-on wave propagating with velocity c from a 
compact object [120]. Inside the switching-on wave, the 
flow rapidly becomes stationary, in correspondence with 
the analytical solution. 



poloidal velocity of particles matches the fast magne- 
tosonic wave velocity goes formally to infinity (and the 
inner surface, oppositely, coincides with the black hole 
horizon). But we have seen that it is the conditions on 
the fast magnetosonic surface that fix the longitudinal 
current circulating in the magnetosphere. In addition, 
the very necessity of establishing the boundary condi- 
tion on the black hole horizon, i.e., in the causally dis- 
connected region, shows that the physical interpretation 
given above does not relate to the reality [132]. 

To make this point more clear and, in particular, to 
understand the boundedness of stationary solutions (and 
hence the analytical method itself) in studies of the cur- 
rent closure, it is useful to consider the results obtained 
in paper [120] in which the problem was formulated as 
follows. There is a magnetized ball at rest which at the 
moment t = starts rotating with angular velocity Q. 
As a result, the switching-on wave starts propagating 
from the ball with velocity c, so that the magnetic field 
remains unperturbed beyond it and electric currents are 
absent, while inside the switching-on wave (and this is 
a very important result) the solution rapidly approaches 
the stationary transonic regime, which is in full agree- 
ment with the analytical solution. Thus, the assumption 
of the stationary solution for longitudal currents flowing 
actually along magnetic surfaces is confirmed. 

As for the current closure, no current closure as such 
happens at all in the ideal case where the outfiow oc- 
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curs in a vacuum. This is related to the fact that in the 
switching-on wave the flow is time-dependent (Fig. 11), 
so there divj 7^ (S S Komissarov, private communica- 
tion). In reality, the current closure will take place on a 
shock wave which must necessarily emerge in the region 
where the supersonic switching-on wave collides with the 
surrounding medium. In any case, however, the ambient 
medium for transonic fiows cannot infiuence the magni- 
tude of the longitudinal current for r < rp and, hence, 
affect the central engine energy release. As soon as the 
switching-on wave crosses the singular surfaces (they are 
also shown in the left lower corner of Fig. 11), the longi- 
tudinal current fiowing in the magnetosphere stops de- 
pending on time. For this reason, one can indeed con- 
sider in the framework of the stationary approximation 
that the electric current closure occurs at infinity, as is 
usually assumed. Thus, transonic flows are signiflcantly 
different from subsonic ones, when the electric current 
circulating in the magnetosphere is determined by the 
conductivity of the boundary of the region occupied by 
plasma (see, for example, Ref. [133]). 

Accounting for the nonzero mass of particles allows us 
to clarify the situation with the 'boundary condition on 
the horizon', and thus with the mechanism of energy re- 
lease by a black hole. Indeed, as shown in Fig. 10, the 
fast magnetosonic surface for nonzero masses of parti- 
cles is located above the black hole horizon. Therefore, 
the critical condition should also be set here, which is 
definitely located in the region casually connected to the 
outer space. The black hole horizon will be located in 
the region of the supersonic fiow and, hence, cannot af- 
fect the properties of the fiow. As a result, the additional 
critical condition must also be kept in the force-free limit 
Trip — >■ 0, when the fast magnetosonic surface, as noted 
above, formally coincides with the event horizon. It is 
not then surprising that this limit of the critical condi- 
tion on the fast magnetosonic surface exactly coincides 
with condition (74) [10]. Thus, the boundary condition 
on the horizon (74), which was necessarily used in the 
force-free approximation, represents a relict of the criti- 
cal condition on the fast magnetosonic surface. 

Correspondingly, it also becomes clear how to inter- 
pret the infinite time retardation near the black hole 
horizon. Indeed, the time it takes for the plasma to 
reach the horizon must be infinite from the point of view 
of the remote observer. As the time of existence of the 
black hole (the surrounding plasma, etc.) is finite, the 
remote observer will register a 'switching-on wave' corre- 
sponding to the very initial stages of the existence of the 
central engine over the black hole surface. This was, in 
fact, one more argument in support of the fact that it is 
incorrect to set any boundary condition on the black hole 
horizon [132]. However, as seen by the example of the 
outflow (see Fig. 11), to form the current system that 
fully determines the central engine power it is sufficient 
to wait a finite time until the plasma intersects the fast 
magnetosonic surface. This also implies, inter alia, that 
after a finite time the switching-on wave will turn out to 
be in the supersonic region of the flow and, hence, can- 
not affect the magnetosphere structure. This conclusion 
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E=-V ^B/c 




Figure 12: The appearance of the electromotive force in a 
circuit at rest relative to a rotating black hole immersed 
in an external magnetic field. The electric field directions 
are shown for observers at rest. 



was numerically confirmed by Komissarov [134]. 

If this is indeed the case, however, another source of 
EMF should be found, since the black hole itself now 
can no longer be the source of extraneous forces and, 
consequently, serve as a battery. It turned out that the 
appearance of EMF in the black hole magnetosphere is 
related to the Lense-Thirring (the frame-dragging) ef- 
fect, which appears due to rotation of the black hole. 
Indeed, according to general relativity, the space itself 
in the vicinity of a Kerr black hole starts rotating with 
angular velocity uj (AA). Only in the reference frame 
rotating with this angular velocity will an observer not 
register a precession of gyroscopes. On the contrary, 
noninertial forces will appear in the laboratory frame at 
rest with respect to remote observers and they can be 
detected. It is the frame-dragging effect that leads to 
the appearance of the electromotive force in the black 
hole magnetosphere. 

Indeed, as in the case of any body moving in the mag- 
netic field, the 'motion of space' relative to the observer 
at rest will produce the electric field E = —V x B/c, 
where now V = w x r is the velocity of the body at 
rest relative to the preferred reference frame. Here, it is 
important that the angular velocity uj be different at var- 
ious distances from the black hole. Therefore, as shown 
in Fig. 12, the circulation of the electric field in a cir- 
cuit will be nonzero even if the electric circuit is at rest 
relative to the black hole, i.e., when the magnetic fiux 
through the contour remains constant (of course, this 
state is possible only above the ergosphere). It is the 
motion of space through the circuit that generates the 
electromotive force. We see here that the 'electric bat- 
tery' will be located above the black hole horizon, and 
possibly outside the ergosphere. 

Finally, as already stressed, for positive energy fiux 
from the central engine the electric currents in the black 
hole magnetosphere can flow only in the case where the 
plasma generation mechanism is operating above the 
horizon. Electron-positron pair production in the col- 



lision of hard gamma-ray quanta emitted from the ac- 
cretion disc surface was already discussed above. Here, 
we should note that, as in the Penrose effect, one par- 
ticle should fall into the black hole, and another par- 
ticle should escape to infinity. It should be recalled 
that the Penrose effect has its origins in the remark- 
able property of rotating black holes — the relativis- 
tic mass defect can exceed 100% inside the ergosphere 
rg < r < rc = M 4- y/AI'^ - cos^ 6* (G = c = 1) and, 
hence, above all the horizon surface [14]. Therefore, it 
becomes evident that the Blandford-Znajek effect is in 
fact the electromagnetic realization of the Penrose pro- 
cess. The difference is that it concerns not charged par- 
ticles themselves but the electromagnetic field they in- 
duce. In other words, the spin-down of the black hole is 
not due to electric currents flowing over the horizon but 
due to the negative electromagnetic energy fiux falling 
onto the black hoi. It is this role of the ergosphere inside 
which the relativistic energy of any material bodies (in- 
cluding the electromagnetic field) can become negative 
that is great. Such interpretation seems now to be the 
most likely, and most researchers involved in these stud- 
ies tend to accept this interpretation (see, for example, 
book [135]). 

Let us summarize. We have shown that the critical 
conditions on the fast magnetosonic surface mostly de- 
termine the energy release of the central engine. This 
surface serves as a valve that determines the magnitude 
of the longitudinal current circulating in the magneto- 
sphere. Like the usual sonic surface in hydrodynamics, 
it separates subsonic and supersonic parts of the flow. 
As a result, the longitudinal current for transonic flows 
must be determined not by the external conditions but 
by the condition of smoothly crossing the singular sur- 
faces. Even more complicated is the situation in a black 
hole magnetosphere, where both the current and the an- 
gular velocity fip should be deduced from the critical 
conditions on the singular surfaces. In consequence of 
this, it is these critical conditions on the singular sur- 
faces that will determine the central engine power. 



4 Theoretical predictions 
4.1 Collimation 

4.1.1 The force balance across the flow. Let now 

analyze the results of the analytical theory. First of 
all, we should try to formulate some general proper- 
ties of a magnetized outflow, which must show their 
worth at large distances from the central object. As 
already noted, we are primarily interested in transonic 
flows, in which the flow at large distances is supersonic. 
Moreover, we shall unconditionally assume that the so- 
lution can be continued to inflnity. This is possible, as 
we have seen, only if the longitudinal electric current 
is sufficiently large. As we have already stressed, ther- 
mal effects can almost always be neglected. As a conse- 
quence, Bernoulli equation becomes a fourth-order alge- 
braic equation with respect io , which in many cases 
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Figure 13: The structure of magnetic surfaces for a non- 
relativistic plasma outflow in the split monopole mag- 
netic field [120]. The decoUimation in the region of the 
reverse bulk current near the central engine is clearly 
seen. However, the redistribution of longitudinal cur- 
rents occurs at far distances, so that the current fiowing 
out is concentrated near the rotational axis, and the re- 
verse current flows near the equatorial plane. 



unlike particles, can intersect the magnetic surfaces. In- 
deed, it is easy to check that, for an almost radial flow 
at large distances, the right-hand side of equation (75) 
should decrease as r~^, while the numerator on the left- 
hand side decreases as r~^. As a result, for equality (75) 
to hold, the radius of curvature of the magnetic surfaces, 
Rc, should increase as r. But such a behavior cannot be 
realized at mathematical infinity [137], so for r — > cx) 
only the right-hand side of equation (75) is leading. 

In the long run, we come at first glance across the para- 
doxical result that the current density in the magneto- 
sphere must vanish at large distances [138]: 



hi = 0. 



(76) 



In fact, this simply means that at large distances al- 
most all outflowing longitudinal current must be con- 
centrated near the rotational axis. Indeed, as shown in 
Ref. [139], a cylindrical region containing almost all out- 
flowing current must inevitably appear near the axis|3 
This behavior was later verified by numerical simula- 
tions [120]. We shall consider this point in more detail 
in Section 4.1.3. 

In the relativistic case, the electric force p^E and the 
component of the Ampere force related to the longitu- 
dinal current j|| x will mostly act on the outfiowing 
plasma. However, as follows from equation (60), these 
forces almost mutually balance each other. Therefore, in 
addition to the bulk forcc0 



allows us to write out rather simple analytical asymp- 
totic solutions. 

For simplicity, let us start from the nonrelativistic 
case. By analyzing the leading terms in the Grad- 
Shafranov equation (41), it is possible to show that the 
force balance equation can be written in the form [136] 



1 



-J\\B. 



(75) 



where Rc is the radius of the magnetic field line in the 
poloidal plane. In other words, the equilibrium must 
be established due to the balance between the centrifu- 
gal force pv^^/Rc and the Ampere force related to the 
longitudinal electric current j|| . It is evident that the 
Ampere force, depending on the sign of the longitudinal 
current, can both collimate and decollimate the flow. In 
this case, the collimation must occur close to the jet axis, 
while the decoUimation occurs at its periphery. Exactly 
this behavior was obtained both analytically [111] and 
numerically [120] (Fig. 13). It should be stressed that 
here we do not consider how strongly the magnetic sur- 
faces can be bent and only investigate their form. 

A similar picture, however, cannot be realized up to 
very large distances from the central engine. 'Natura 
abhorret vacuum', and the diverging magnetic surfaces 
inevitably must start collimating. It turned out that this 
is possible exactly because the longitudinal current / is 
not the integral of motion, and hence the electric current. 



jpol 




(77) 



it is necessary to take into account the bulk centrifugal 
force 

nm(?^ -f S/c 



Rc 



(78) 



[S « cB-^/At: is the Foynting vector) and the Ampere 
force J^jtor ~ ji^ X Bp/c pertaining to toroidal current. 
As a consequence, the Grad-Shafranov equation in the 
limit r ^ r-p can be conveniently rewritten in the form 
[136, 138, 140] 
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(n ■ 



.) = 0, (79) 



where n = V5'/|V^'|. Notice that the Foynting vector 
contributes, in addition to the contribution from parti- 
cles, to the centrifugal force. This is related to the fact 



* Strictly speaking, this terminology corresponds to the case of 
f2B < 0, where pQj > 0. For the opposite orientation, the current 
near the rotational axis will flow towards the central engine. 

^The corresponding formula (4.227) from monograph [9] has 
the incorrect sign 
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that, as noted above, both particles and the electromag- 
netic energy propagate along the magnetic surfaces. 

4.1.2 The collimation mechanism. Thus, the form 
of the magnetic surfaces close to the rotational axis de- 
pends on the balance between the collimating Ampere 
force arising from the longitudinal electric current (paral- 
lel currents are attracted) and the decollimating Ampere 
force related to toroidal currents. Thus, the question as 
to whether the collimation will be effective depends on 
the magnitude of the longitudinal current. A series of ex- 
act solutions [110, 111], which were possible to obtain by 
analyzing small deviations from the monopole magnetic 
field, showed that for nonrelativistic jets the collimation 
is large even close to the fast magnetosonic surface. This 
property is also confirmed by numerical modeling [120, 
141] (see also Fig. 13). 

In the relativistic case, where, as we remember, the 
longitudinal current is close to the Goldreich current 
(zo ~ l)j an almost full compensation of these two forces 
takes place. In particular, the balance is met exactly in 
the force-free approximation and in the split monopole 
magnetic field [142], so the vacuum monopole solution 
remains exact up to infinity for the magnetosphere filled 
with the plasma as well (Fig. 14). The current I here 

takes the form I{9) — I^' sin^ 9, where 



which exactly coincides with the Goldreich current 
(j|| = PGjc). 

For massive particles (and again for the split monopole 
magnetic field), the longitudinal current determined 
from the condition of smoothly crossing the fast mag- 
netosonic surface will differ from the Goldreich current 
only by a factor of the order a~^^^ [143]. As a result, 
the perturbation of the magnetic flux function (54'/\l/ in 
the asymptotically remote region r ^ rp will increase 
logarithmically slowly [143,144]: 

In other words, the current turns out to be only insignif- 
icantly larger than the critical one, which leads to a van- 
ishingly small collimation. Correspondingly, the particle 
energy also increases very slowly: 

7«ai/3lni/3(^-^y (82) 

We stress that above we have considered the proper 
collimation, i.e., that due to bulk currents. However, 
the collimation, generally speaking, can be produced in 
the source itself. In Fig. 15, the flow obtained as a small 
perturbation of the force-free solution is shown, but for a 
parabolic field [145]. Such a field may be generated in the 
accretion disc as well [5]. For not too small {9 ^ 7in/'^) 
and not too large (6' <^ a^^/'^) angles, the location of the 




Figure 14: The force-free monopole solution found by 
Michel [142], in which the electric field Eg exactly equals 
the toroidal magnetic field . The contour arrows show 
the direction of the poloidal current. 




Figure 15: The location of the Alfvenic (A) and fast 
magnetosonic (F) surfaces in a parabolic magnetic field. 
Such a field can also be generated in the accretion disc 
[145]. 
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fast magnetosonic surface rp and the value of 7(rF) are 
given by the expressions 



(83) 
(84) 



respectively, which again correspond to the estimates of 
the particle energy, 7 ~ a^^^, and the axial distance, 
rp sin 9 ^ a^^^Ri^ obtained above (see Ref. [145] for 
more detail). But in this case the account for nonzero 
particle masses also only slightly perturbs the force-free 
solution. On the fast magnetosonic surface, for example, 
one finds 

r—r-p 



Thus, we can conclude that the proper coUimation is 
impossible in the relativistic case. That is, the bulk 
collimation due to longitudinal currents flowing in the 
magnetosphere is impossible. Therefore, the magnetic 
surfaces can be collimated either by a special choice of 
currents in the source itself, or due to the interaction 
of the supersonic wind with the surrounding medium. 
Clearly, a sufficiently extended dense disc is required in 
the first case. Because of this, for radio pulsars where 
there is definitely no such disc, one usually assumes that 
the flow at distances r ^ Ri, must be radial. 

Unfortunately, it is not clear at present which mech- 
anism is actually responsible for the collimation of rela- 
tivistic jets. However, if one assumes that the collima- 
tion is indeed caused by the presence of the surrounding 
medium, it becomes possible to estimate the transverse 
size Tjet of the jets. Indeed, by assuming that the pres- 
sure of the poloidal magnetic field in the jet is close to 
the ambient pressure Pext (this estimate is valid by as- 
suming the total electric current in the jet to be zero), 
we obtain from the condition of the magnetic flux con- 
servation that 



net 



R 



Bl 



1/4 



(86) 



where R and B-^ are the radius and the magnetic field of 
the compact object, respectively. For example, for active 
galactic nuclei (Bin ~ 10* G, i? ~ 10^^ cm) we have 



Ipc, 



(87) 



which exactly corresponds to the observed transverse 
sizes of jets [121. Correspondingly, for young stellar ob- 
jects (Bin ~ 10^ G, i? ~ 10" cm) we have rj ~ 10^^ cm, 
which again is in agreement with observations. Thus, the 
external medium apparently must significantly affect the 
collimation of jets. 

4.1.3 The dense core. As we have seen, near the 
axis the bulk force J^jpoi relatc;(l to the poloidal current 
is always directed toward the rotational axis, and in the 
region of the current closure, away from the axis. Then, 



for example, close to the cylindrical core where the cur- 
vature of the magnetic surfaces is small and, hence, the 
centrifugal force J-"cont (78) can be neglected, the equilib- 
rium can be established only if the force J^jtor is directed 
away from the rotational axis. But the poloidal mag- 
netic field in this case must decrease with axial distance. 
Then, the density of the outflowing plasma will decrease 
along with the poloidal field as well. Exactly such a 
behavior of the solution was illustrated in Fig. 13. 

It turned out that the dense core for the nonrelativistic 
supersonic flow will exist for both strongly and weakly 
magnetized flows close to the central engine, irrespective 
of the ambient pressure [146]. The radius of the core will 
correspond to such a distance from the axis, at which 
the toroidal field matches the poloidal one. A straight- 
forward calculation shows that in both cases the longi- 
tudinal current in the core region will be j = (c/uin)jGJ, 
i.e., io = c/viri, and so one arrives at 



^core — 



The magnetic flux in such a core must be a significant 
fraction of the total magnetic flux. Such a configura- 
tion is none other than the z-pinch well-known in plasma 
physics [147] (the question of stability will be briefly dis- 
cussed below). 

As for the relativistic flow, the appearance of the dense 
core can be balanced here by the electric force which, 
as we have seen, significantly weakens the force J^jpoi- As 
a result, the answer is significantly dependent on the am- 
bient pressure [9]. If the relativistic jet is surrounded by 
a medium with the total gas and magnetic prc^ssure Pext 
above some limiting value Pmin, the dense core does not 
form at all. In this case, the poloidal magnetic field in 
the jet will be approximately constant: Bp/Sn ~ Pext- 
For convenience sake, we shall express below the lim- 
iting value Pext through the equivalent magnetic field 
(P^i„/87r = Pmin). Here, one obtains the relationship 



B„ 



-BiRi.), 



(89) 



where B{Rl) = fl^^tot/Trc^ is the characteristic mag- 
netic field on the light cylinder. Correspondingly, the 
density of the outflowing plasma will be constant, too. 

If the ambient pressure is sufficiently low, so that 
-Pext < Pmin, then, as in the nonrelativistic case, the 
dense core is formed in the center of the outflow; the 
radius of the core must exceed that of the light cylinder: 



^ 7inPL- 



(90) 



It is easy to check that on the core boundary the energy 
flux density of the electromagnetic fleld matches the par- 
ticle energy density. The appearance of such a cylindrical 
jet was predicted already many years ago in many papers 
[138, 148, 149], but the magnetic flux in this core was 
determined only quite recently [145,146]. The magnetic 
fleld near the axis was found to only very weakly (log- 
arithmically) depend on the ambient pressure, so with 
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Figure 16: Longitudinal magnetic field Bz (G) as a func- 
tion of the axial distance x — Q,uj/c for a — 10^ and 
different Mach numbers Ai on the jet axis, which was 
obtained as a solution of the one-dimensional problem 
[146]. The axial magnetic field only slightly deviates 
from -Bmin when the ambient pressure changes by sev- 
eral orders of magnitude. 



good accuracy one can set Score = ^min- As a result 
the magnetic flux = '"'''corc^min within the cor( 

turns out to be much smaller than the total flux: 



^«1. 



(91) 



Here, the poloidal magnetic field and the density of mat- 
ter for Xr = tn/rcorc > 1 must behave as power functions 



(lab) -fe2 



(92) 
(93) 

respectively. As the ambient pressure decreases, the ex- 
ponents ki and k2 gradually increase; however, their dif- 
ference remains approximately constant, viz. 



ki - k2 



0. 



(94) 



As a result, if we have fci w ^2 ~ for -Pcxt > ^min 
(i.e., the poloidal magnetic field and the outfiowing 
plasma density are constant in the cross section), for 
ambient pressures corresponding to the magnetic field 

^ext 



Bcq, where 



1 



(95) 



we have, in contrast, A:i w ^2 ~ 1. Under such ambient 
pressures, as we shall see, the contribution from particles 
to the energy flux becomes dominant in the entire volume 
of the jet. 

In Fig. 16, the behavior of the poloidal magnetic field 
is shown in the double-logarithmic scale as a function of 
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Figure 17: Longitudinal magnetic field profile as a func- 
tion of the axial distance ^ oc tu for different values of 
the ambient pressure, which was obtained as a solution of 
the two-dimensional time relaxation problem [126]. The 
power-law dependence corresponds to the analytical es- 
timate (92). 



the distance to the axis [146] (see also Ref. [150]). It 
was obtained as the solution to the system of two ordi- 
nary differential equations (^.14) and {A.27) to which, 
as shown in the Appendix, the Grad-Shafranov equa- 
tion can be conveniently reduced in the cylindrical case. 
The Mach number on the rotational axis served as the 
parameter for different curves. It is evident that the 
poloidal magnetic field at distances exceeding Tcorc in- 
deed starts decreasing as a power law. At the boundary 
of the jet it changes by several orders of magnitude, while 
the magnetic field on the axis changes only several-fold 
times. This behavior was recently confirmed by numeri- 
cal two-dimensional calculations, too [126] (Fig. 17). It 
is important here that in the last case the very formu- 
lation of the problem was principally different (the time 
relaxation problem was solved rather than a stationary 
problem) . 

In relativistic jets there can be one more regime which 
is impossible in nonrelativistic jets [151]. In Fig. 15, 
it corresponds to cross sections located sufficiently close 
to the equator, when the flow in the jet center must be 
subsonic. By rewriting the condition z < {a/jin)R'L in 
terms of an ambient pressure, we obtain the inequality 



Pext > 



Btt, where 



Br 



= ^S(i?L) 



(96) 



In this subsonic flow region must inevitably be 

formed in the inner parts of the jet with zu < rg, where 



-,1/2 
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(97) 
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At last, the subsonic flow is established in the entire 
jet volume for higher ambient pressures corresponding 
to the magnetic fleld pressure on the fast magnetosonic 
surface. 

4.1.4 The current closure. As wc have understood, 
according to the most wide-spread point of view, it is the 
electric current flowing along collimated outflows, which 
is responsible for the main jet energy release. Here, the 
following theorem can be formulated. 

A stationary cylindrical jet with finite magnetic flux 
\l/tot can be formed either with the nonzero full electric 
current /(^ftot) ^ or in the presence of a surrounding 
medium with nonzero pressure. 

At first glance, these two variants fully contradict each 
other. However, this is not actually the case. As wc have 
seen, a cylindrical core (which contains a significant part 
of the outflowing current) must be formed near the jet 
axis both in the relativistic and nonrelativistic cases. If 
the compact object were solitary in the Universe, the 
reverse current would indeed return near the equatorial 
plane. The closing of the current itself would occur in the 
switching-on wave. Thus, if one ignores the current clo- 
sure region and studies only the structure of the central 
region, it is indeed possible to assume that the current 
closure takes place at inflnity. 

In fact, when a jet is immersed into a medium with 
a finite pressure, the reverse current, as one usually as- 
sumes, must flow along the boundary of the cocoon that 
forms due to the interaction of the supersonic wind with 
the external medium. The boundaries of the correspond- 
ing cocoons can be well seen in Figs 1 and 4. It is the 
force balance at the cocoon boundary that allowed the 
magnetic fleld Smin close to the jet axis to be determined 
[145]. 

4.1.5 The stability. In conclusion, it is necessary to 
briefly consider the problem of the stability of jets. The 
nonrelativistic z-pinches observed in the laboratory are 
known to be strongly unstable with respect to constric- 
tions and screw modes [147, 152,153]. Therefore, the 
problem of jet stability has been widely discussed in the 
literature [154-160]. As shown in Fig. 4, however, jets 
from young stars at large distances indeed are similar to 
a sequence of flying blobs rather than to a regular flow. 
Only at small distances can the flow be considered quite 
regular (see Fig. 5). 

In recent years, special laboratory experiments have 
been carried out under conditions maximally similar to 
those in the nonrelativistic jets [161, 162]. In partic- 
ular, the plasma velocity was as high as 100-400 km 
s~^ (which is comparable to the plasma flow velocity in 
jets from young stars), and the total current was about 
1 MA. In these experiments, a strong instability leading 
to rapid fragmentation of the flow into individual blobs 
was also observed. 

Nevertheless, it is not obvious that the instability of 
laboratory pinches can be considered as an undisputable 
evidence for the immanent strong instability of astro- 
physical jets. The point is that astrophysical jets are 
always 'specially prepared', since they come out from 



the quasispherical subsonic flow region. As a result, the 
plasma density and the longitudinal magnetic field pro- 
files near the jet base turn out to be close to the equilib- 
rium ones. In laboratory experiments, in contrast, the 
initial plasma density is usually very different from the 
equilibrium value [153]. 

As for relativistic jets (which at large distances also 
frequently show an irregular structure), they proved to 
be more stable [163 165]. The recent numerical simu- 
lations [166] (where the jet was found to be stable af- 
ter more than 1000 rotations of the central engine) con- 
firmed this conclusion. Thus, there are no doubts now 
that the nonrelativistic Kruskal-Shafranov stability cri- 
terion [153] 



where L is the length of the jet, cannot be applied to 

relativistic flows. Unfortunately, the limits of the present 
review do not allow us to discuss this most important 
point in more detail. 

4.2 Acceleration 

4.2.1 The acceleration mechanism. First of all, let 
us consider the acceleration mechanism itself. It is conve- 
nient to start from expression (43) for the toroidal veloc- 
ity of the non- relativistic flow. It is evident that in the 
subsonic region Ai^ < 1 the flow velocity corresponds to 
the precise corotation: 

« Qf-cu. (99) 

In other words, particles can be considered as beads on 
a wire that determines their angular velocity of rotation. 
This situation is quite understandable because, as can be 
easily checked, the energy density of the magnetic field 
within the Alfvenic surface exceeds the plasma energy 
density, so it is the magnetic field that controls the mo- 
tion of particles. In fact, the magnetic field plays the 
role of a slingshot that provides the constant angular ve- 
locity of the plasma rotation. Therefore, the velocity of 
particles linearly increases with increasing axial distance. 
Then, the following theorem can be formulated. 

In the nonrelativistic limit, the smooth crossing of the 
fast m,agnetosonic surface is possible only if the particle 
energy flux on this surface is at least one-third of the 
total energy losses. In other words, the transonic non- 
relativistic flow in the supersonic region must already be 
effectively accelerated (q ^ I). 

Indeed, we can set _En « ilF//27rcr^n for strongly mag- 
netized flows near the surface of the central engine, which 
immediately yields « 0|rp for / « zq/gj- Moreover, 
it is easy to show that the poloidal velocity near the 
singular surfaces also becomes on the order of Ofi^f, too. 

By this means the particle acceleration mechanism in 
the strongly magnetized wind is similar to that by which 
a slingshot accelerates a stone. This becomes possible 
exactly due to the dominant effect of the poloidal mag- 
netic field. However, this acceleration stops at distances 
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fp ~ ^'Ai since there the magnetic field energy density 
becomes smaller than the plasma particle energy den- 
sity. Additionally, the toroidal magnetic field beyond the 
Alfvcnic surface becomes stronger than the poloidal one, 
so particles start sliding relative to the magnetic field 
lines. As a result, the flow passes to another asymptotic 
behaviour: w Qpva\/w.. The plasma kinetic energy 
in this region will be mainly due to their poloidal velocity 
component. 

As for the particle acceleration in the ultrarelativistic 
limit, near the fast magnetosonic surface, as we have 
seen, the particle energy flux must be much smaller than 
the electromagnetic energy flux. So the question arises as 
to whether it is possible to effectively accelerate particles 
beyond the sonic surface. Surprisingly, the force balance 
across the flow should be considered once again to answer 
this question. 

4.2.2 EfRciency. So, let us come back to equa- 
tion (79). We have already mentioned that the particle 
energy will be completely determined by the drift mo- 
tion beyond the light cylinder w ^ Ri,. Therefore, using 
Eqn (59) we can write out the expression for the Lorentz 
factor of particles in the form 
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(100) 



lEl 



r-Bp, where 



Now, making use of the relation 

again Xr — ^Iw/c, we immediately come to the conclu- 
sion that when the second term in expression (100) can 
be neglected, the particle Lorentz factor must approach 
the following asymptotic solution: 



(101) 



which, as we shall see, is universal enough. If the cur- 
vature of the magnetic field lines is appreciable, then, in 
contrast, we can neglect in formula (79) the second term 
corresponding to the bulk force j^Bp/c. As a result, by 
comparing the corresponding terms in the force balance 
equation (79) for the strongly magnetized (i.e., Poynt- 
ing dominated) flow, we arrive at another asymptotic 
solution [167]: 



(102) 



where C ~ 1. Moreover, making use of Eqn (59) and the 
equilibrium condition (79), we can write in the general 
case the relationship [109] 
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2 C^Rc' 



(103) 



Here, the value of C can be exactly determined for 
strongly eollimated flows and Op = const [109,150]: 



Vs. 



(104) 



Simply speaking, the Lorentz factor will be determined 
by the least of the two values giving by expressions (101) 
and (102). 
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Figure 18: The growth of a particle's Lorentz factor 7 
with distance from the equatorial plane z [168]. The 
right curve corresponds to strong collimation (A; > 2). 
For weaker collimation (left curves), the particle accel- 
eration at large distances becomes less eff'ective. 



Thus, the choice between asymptotic solutions (101) 
and (102) must be determined by how bent the mag- 
netic surfaces are. It is easy to show that the parabolic 
magnetic field in which the field line at a large dis- 
tance from the central source is given by the equa- 
tion z(zu) cx zu"^ corresponds to the terminating case 
[126,146]. Indeed, as the curvature radius can be defined 
as i?c = [{z')'^ + l]^/'^/z", in the limit z' = dz/dw > 1 
for the magnetic surfaces specified by the relationship 
z{w) (X zu'', Eqn (102) gives the energy of particles mov- 
ing along the magnetic field line: 



7 a ro' 



fe-i 



(105) 



where w is, in this case, the current axial distance of a 
particle. At fc = 2, the acceleration efficiency determined 
from expressions (101) and (102) is the same. Thus, if 
the magnetic surfaces are eollimated more strongly than 
those for the parabolic field (i.e., > 2), the curvature 
of the magnetic surfaces can be neglected at large; dis- 
tances, and the energy of the particles will be determined 
by expression (101). If the flow is poorly eollimated (i.e., 
1 < fc < 2), the particle acceleration will be less effec- 
tive and one should use expression (102). The numerical 
simulations [124, 168] fully confirm the picture presented 
here. As shown in Fig. 18, the acceleration of a particle 
moving along the magnetic field line indeed follows the 
law 7(2;) oc z^/*^ for strongly eollimated flows, in corre- 
spondence with asymptotic behaviour (101), while the 
acceleration for poorly eollimated outflows at large dis- 
tances, in full agreement with expression (105), becomes 
less efl[ective. 
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We thus arrive at the most important conclusion that 
the acceleration efficiency of particles in a supersonic ul- 
trarelativistic wind is determined by the degree of col- 
limation of the magnetic surfaces. The plasma can be 
effectively accelerated only in the case where the mag- 
netic surfaces are collimated more strongly than those for 
the parabolic field. In this event, one can use asymptotic 
solution (101) which shows that the acceleration mech- 
anism, in fact, again is due to the slingshot effect (the 
larger the axial distance, the higher the energy). But 
in this case, too, the total transformation of the elec- 
tromagnetic energy into the particle flux energy, 7 « ct, 
can take place only if the jet transverse size exceeds the 
value 

rcff = <yRL- (106) 

In particular, the fraction of energy carried by particles 
for Pcxt < B^q/^T^ can be determined from the simple 
relationship [10,151] 



part 



1/4 



(107) 



If the coUimation is poor, the acceleration will be inef- 
fective, since particles begin sliding from the magnetic 
field lines. In this case, a much higher transverse size of 
the jet is required for accelerating particles to limiting 
energies: 

reff = ai/('=-i)i?L. (108) 

This dependence was also numerically confirmed [169]. 
In particular, the acceleration becomes actually impossi- 
ble for the split monopole magnetic field k = 1 (which, as 
we see, is a special case that requires separate considera- 
tion), since, as follows from estimate (82), the condition 
7 ft! cr is reached only at exponentially large distances 
from the compact object. On the other hand, it easy to 
check that for fc > 1 both for effective and for ineffective 
acceleration mechanisms the condition 71? ~ 1 will be 
satisfied, where ~ ro/z is the characteristic opening 
angle of the acceleration region. 

Here, it is important to emphasize that the self- 
consistent analysis, in which the magnetic surfaces for 
the split monopole field were not assumed to be exactly 
conical, allowed the determination of the correct struc- 
ture of the flow [143]. For example, the expression (82) 
for the particle's Lorentz factor discussed above is deter- 
mined precisely by the small curvature of the magnetic 
surfaces, and so exactly corresponds to expression (102). 
It should be remembered that for a long time the flow 
in the split monopole magnetic field had been, in fact, 
the only example in which analytical results could be ob- 
tained. Therefore, a strong opinion was that the effective 
particle acceleration beyond the fast magnetosonic sur- 
face is completely impossible. As wee see, this conclusion 
proved to be incorrect in the general case. 

4.3 Subsonic flows 

As discussed above, all observed jets must be transonic. 
It is this property that allowed us to find the longitu- 




Figure 19: The structure of the axisymmetric radio- 
pulsar magnetosphere in the model [98]. The fiow is 
assumed to be radial at large distances. 



dinal current flowing in the magnetosphere and, hence, 
to estimate the central engine power. However, until re- 
cently many models considered subsonic flows, and the 
magnitude of the current was determined from other con- 
siderations [170, 171]. We shall briefly discuss below two 
the most known examples of such subsonic flows. 

4.3.1 Pulsar magnetospheres. Starting in the 
early 1970s, pulsar magnetospheres have been discussed 
mostly in the force-free approximation [107, 142, 172- 
174]. This was based on the fact that the plasma fill- 
ing the neutron star magnetosphere is secondary with 
respect to the magnetic field, and so (at least inside 
the light cylinder) the particle energy density can be 
neglected. The Grad-Shafranov equation (^.28) in the 
force-free approximation (which in this case is simply 
called the pulsar equation) becomes elliptical. There- 
fore, for numerical modeling of the axisymmetric mag- 
netospheres one has to impose an additional condition 
at the external boundary of the integration region [98, 
100 104]. Usually, one chooses the condition of radial- 
ity of the magnetic field lines (Fig. 19). In this case 
it is this additional condition that fixes the longitudinal 
current in the magnetosphere. Thus, it is not surprising 
that this current turns out to be close to the longitudi- 
nal current /m(^) (80) obtained by F C Michel for the 
monopole solution shown in Fig. 14. If the absence of 
the reverse current fiowing along the equator (so that 
the current closing was done by bulk currents only) was 
chosen as the additional condition, the magnetic field 
structure beyond the light cylinder turned out to be sig- 
nificantly different [175]. 

Interestingly, a similar structure also appears in dif- 
ferent models of the inclined rotator. Here, first of 
all, we should highlight the model of the 'rotating split 
monopole' (Fig. 20). In the force-free approximation 
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Figure 20: The magnetic field structure in the 'rotating 
spht monopole' model [176]. In the force-free case, parti- 
cles moving radially at the speed of light can provide the 
formation of the current sheet (the wavy curve) which 
separates magnetic fluxes in the equatorial region. 



(when massless particles move with the speed of light), 
the monopole magnetic field, namely 

^{r,9,ip,t) = ^-0(1 - cos 6*), 

61 < 7r/2 - X cos(^ -nt + ilr/c), (109) 
^{r,0,ip,t) = ^-0(1 -I-C0S6'), 

9>tt/2~x cos((p ~nt + nr/c), (110) 

was found to be also the solution of the problem [176]. 
In this case, the electromagnetic fields inside cones with 
angles 9 < tt/2 — x and n — 9 < Tr/2 — x near the rota- 
tional axis are independent of time and the angle ip and 
coincide with the fields of an axisymmetric rotator, while 
in the equatorial plane the field directions reverse sign at 
the moment when the current sheet intersects the given 
point. A similar structure beyond the light cylinder was 
obtained for the rotating dipole, too [177]. True enough, 
in contrast to the 'rotating magnetic monopole', the to- 
tal longitudinal current here depends on the inclination 
angle Xj although not very strongly. 

The last property can be easily explained. As we have 
seen, the very existence of the MHD wind far away from 
the light cylinder is possible only if the toroidal magnetic 
field is comparable in magnitude to the electric field. 
But this is possible only if a sufficiently large longitudi- 

(A) 

nal current / /j^ fiows in the magnetosphere. Let us 
keep in mind that this value of the current can also be 
found from the condition of the smooth crossing of the 
singular surfaces in the full MHD version. In none of the 
numerical calculations mentioned above were bounds on 
the longitudinal currents fiowing from the neutron star 
surface set. Thus, it is not surprising that the longitu- 
dinal current obtained from the solution of the problem 
considered turned out to be on the order of /m(6')- As 
a result, the energy losses in model [176] were found to 
be independent of the angle x- Energy losses in model 
[177] were also found to be weakly dependent on the in- 



clination angle: 



1 Bln^R'^ 



sin^ x) 



(111) 



However, the following problem emerges here. All the 
theories of particle generation near the neutron star mag- 
netic poles [178-180] unambiguously suggests that the 
longitudinal current density cannot be higher than that 
of the local Goldreich current which, as follows from def- 
inition (9), depends on the inclination angle x 



JGJ 



cosx- 



(112) 



Indeed, the local Goldreich charge density near the mag- 
netic poles for the orthogonal rotator must be (57i?/c)^/^ 
times smaller than in the axisymmetric magnetosphere. 
Thus, the longitudinal current flowing along open mag- 
netic field lines should be correspondingly smaller (for 
ordinary pulsars this factor can be as high as 10^). In 
the 'rotating split monopole' model considered above, 
this problem does not arise, because at any inclination 
angle in the polar magnetospheric regions the current is 
always the same as in the axisymmetric case. Just this 
current provided the necessary toroidal magnetic fleld. 
For the inclined dipole, in contrast, it is necessary to 
additionally assume that the longitudinal current in the 
polar cap regions can be significantly higher than the 
local Goldreich current (A Spitkovsky, private commu- 
nication). 

There is one more problem related to the decrease in 
the longitudinal current density as x ^ 90° . The point 
is that the current losses for the local Goldreich cur- 
rent must decrease as the inclination angle x increases 
[56,107]: 

W,ot = ^^^^^^ocosx■ (113) 

Here, the coefficient f^, — 1.59-1.96 depends only on 
the inclination angle x- It is necessary to stress that in 
addition to the cos x factor (which is related to the scalar 
product Wtot = — Jl • K, where K is the braking torque), 
the significant dependence of the current losses Wtot on 
the inclination angle is also contained in the quantity zq. 
The matter is that in the definition of the dimensionless 
current io — I / Iqj in expression (113), the denominator 
contains the Goldreich current for the axisymmetric case, 
whereas at nonzero angles x the Goldreich charge density 
itself depends on the angle x near the magnetic poles. 
It is logical to expect that the dimensionless current io 
for the inclined rotator will be bounded from above as 
^^^\x) ^ COSX- As a result, the current losses, in 
contrast to relationship (111), must decrease with angle 
X at least as cos^ x- In particular, for x = 90° (when 
cos^ X must be replaced by its characteristic value within 
the polar cap region cos^ X; obtain [56] 



= ci^^ — — I «A. 



(114) 
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Figure 21: The structure of electric currents flowing near 
the magnetic poles of an orthogonal rotator. Currents 
flowing along the sepratrix (thin arrows) that separates 
the regions of the open and closed magnetic field lines 
are adjusted to bulk currents (contour arrows) in such a 
way that the closing surface current is fully concentrated 
within the polar cap. 



Here, za = 1 for the local Goldreich current, and the 
coefficient cj_ ~ 1 now depends not only on the profile 
of the asymmetric longitudinal current, but also on the 
polar cap shape. 

Usually, when discussing this point, the following 
counter-argument against the decrease in current losses 
with increasing x is invoked. In the expression for the 
braking torque, namely 



K = - 



i y [r X [J, X B]]d5, 



(115) 



The surface current Js indeed must decrease as cosx 
as angle x increases. But then the characteristic dis- 
tance from the axis to the polar cap, in contrast, will 
increase as sin^, so ultimately the losses will be weakly 
dependent on the inclination angle x even for the local 
Goldreich current. 

As follows from a more precise analysis [107], how- 
ever, the real structure of the surface currents in the 
polar cap region was ignored in this obvious, at first 
glance, consideration. Referring to Fig. 21, the closing 
currents in fact must be such that the current averaged 
over the polar cap region is zero. As a result, one needs 
to consider higher-order effects (like, for example, the 
neutron star surface curvature effect), when determin- 
ing the radio pulsar spin-down rate. But if the surface 
current averaged over the polar cap region is indeed zero, 
then, as shown in Fig. 21, a surface current comparable 
in amplitude to the bulk current in the magnetosphere 
must flow along the separatrix separating the regions 
of open and closed magnetic field lines. Remarkably, 



in numerical modeling of the inclined rotator [181], re- 
verse currents flowing along the separatrix were indeed 
discovered. Finally, it worth noting that there is no con- 
tradiction between relations (111) and (113), either. As 
we specially emphasized, the approximation expression 
(111) was obtained in Ref. [177] for a flow in which the 
longitudinal current was larger than the local Goldre- 
ich current, which corresponds to the condition iq > 1 
(lA > 1, respectively). 

In any case, studies of the last decade, in our opinion, 
have at last formulated the problem whose solution will 
provide significant progress in the understanding of the 
structure of radio pulsar magnetospheres. The problem 
is whether the plasma generation region in the neutron 
star magnetosphere can provide a sufficiently large lon- 
gitudinal current which is necessary to launch the MHD 
wind from the inclined rotator. If the necessary current 
can be produced, nothing will prohibit the formation of 
the MHD wind in which the main part of the energy will 
be carried by the electromagnetic field. If the generation 
of a current which is significantly larger than the local 
Goldreich current is impossible, then the toroidal mag- 
netic field near the light cylinder will be smaller than the 
poloidal magnetic field. In that case, a light surface will 
inevitably be formed near the light cylinder, on which 
the current closing and particle acceleration up to en- 
ergies 7 ~ cr will occur [56]. Thus, the problem of the 
effective particle acceleration in pulsar winds, which we 
mentioned in Section 2.4, can be solved. 

Interestingly, the possibility of answering this question 
has apparently emerged a short time ago. This test is 
related to the unusual properties of the radio pulsar PSR 
B1931-h24 [182]. It differs from other pulsars in that it 
stays in the active state for 5-10 days, and then its radio 
emission switches off in less than 10 s, and the source is 
not observed during the next 25-35 days. It is important 
that the spin-down rate in these two states be different: 



f^on = -1.02 X 10 
f2off = -0.68 X 10 "s 
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14 



SO that 
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A similar behavior was later observed in the pulsar PSR 
J1832-h0031 (ton ^ 300 days, toe ~ 700 days), with the 
ratio fion/f^off ~ 1.5 again. 

It is logical to assume that the difference in the spin- 
down rates in these pulsars is simply associated with the 
fact that the spin-down in the switch-on state is due 
to the current losses, and in the switch-off state, when 
the magnetosphere is not filled with plasma, is due to 
magnetodipole radiation [183, 184]. Then, making use 
of equations (3) and (113) we obtain 



off 



^ cot^X, 



(119) 



which yields the reasonable inclination angle x ~ 60°. 
On the other hand, using expression (111) [177] for the 
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switch-on state, we arrive at the relationship 

n^^Sil + su^ (120) 
2 sm X 

Clearly, this quantity cannot be equal to 1.5 for any x- 
If this interpretation of observations holds true, this im- 
plies that the longitudinal current in the magnetosphere 
indeed does not exceed the local Goldreich current. 

It should be emphasized that we have assumed above 
that the magnetospheric plasma fully screens the mag- 
netodipole radiation of the neutron star. This conclu- 
sion, which was formulated for the first time in paper 
[107], seems now to be directly confirmed, since both in 
the model [176] and in the numerical calculations [107] 
there are no alternating electromagnetic fields decaying 
as 1/r. 

4.3.2 'The magnetic tower'. 'The magnetic tower' 
is the model of collimated jets proposed by D Lynden- 
Bell in 1996 [128] (Fig. 22) and later on widely dis- 
cussed in relation to both relativistic and nonrelativistic 
sources [171, 185-187]. It is based on the assumption 
that there exists an intensive wind outflowing perpen- 
dicular to the accretion disc and stretching the magnetic 
held lines along the rotational axis. Here, one usually 
assumes that the initial magnetic held was quasi-dipole, 
i.e., it consisted of magnetic field lines with one end 
frozen in the central star (or in the inner regions of the 
accretion disc) and another end frozen in the outer parts 
of the disc. Such a cylindric flow cannot intersect sin- 
gular surfaces since the magnetic held lines remain at 
a constant distance from the axis of rotation. That is 
why, the longitudinal current in this model will be de- 
termined exactly by the differential rotation that leads 
to the magnetic fleld line twisting. Because of this, as in 
the case shown in Fig. 11, the magnetic tower top will 
propagate upwards by gradually increasing the volume 
occupied by the twisted magnetic fleld, while a station- 
ary configuration restricted by the ambient pressure will 
be formed at smaller distances. 

There are two additional important properties that 
distinguish the magnetic tower model from the transonic 
flows we are considering. First, if the magnetic fleld lines 
are anchored in the accretion disc, a conflguration with 
almost zero total magnetic flux will form during the out- 
flow (see Fig. 22). In other words, the direction of the 
poloidal magnetic field on the periphery of the jet will 
be different from that near the jet axis. Second, when 
the magnetic field loops do not extend beyond the light 
cylinder and, hence, do not open, the energy will con- 
tinue being transferred along the magnetic fields lines. 
But this means that the energy will be carried away from 
the central engine only along the rotational axis, while 
on the jet periphery the energy flux will be directed back 
to the accretion disc [186,187]. 

If such a configuration were stationary, the reverse en- 
ergy flux would be exactly the same as the energy flux 
outgoing from the central engine closely to the rotational 
axis. However, we have seen that the equipotenital con- 
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Figure 22: The magnetic field structure in the model 
of the 'magnetic tower' that can exist under sufficiently 
strong ambient pressure P [128]. The energy is trans- 
ferred from the central engine along the rotational axis, 
and back to the accretion disc along the jet periphery. 



dition is violated in the switching-on wave, so the reverse 
energy flux turns to be smaller than the escaping flux. 

Thus, in the very setup the problem of the magnetic 
tower formation is different from that describing tran- 
sonic flows. The longitudinal current determining the 
energy losses in no way relates to the critical conditions 
on the singular surfaces which, as we specially stressed 
above, should unavoidably appear in all real compact 
sources (see also Section 2). Hence, apparently, the mag- 
netic tower model cannot correspond to the reality. The 
results of numerical simulations also support this con- 
clusion. When the flow intersected the singular surfaces, 
the magnetic field lines became open and the energy flux 
was directed outward from the central engine in the re- 
gions of both outgoing and incoming magnetic fleld lines 
[188]. On the other hand, in the case where the flow re- 
mained subsonic in numerical simulations, the magnetic 
tower formation was indeed observed [50,189]. 

5 Estimation of parameters 

Thus, we have seen that at present it has turned out to 
be possible to understand many key points related to the 
formation and the internal structure of collimated out- 
flows. In doing so, we have managed to flnd several key 
parameters that determine the basic physical properties 
of ejected matter. First and foremost, these include the 
magnetization parameter a, the multiplicity parameter 
A, and the initial velocity v-m (the Lorentz factor 7in) 
of the outflow. We shall try below to understand how 
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precisely these parameters can be estimated for the ob- 
served jets. 

It should be emphasized from the very beginning that 
each time it is necessary to clearly separate which colli- 
mated outflows are being considered. Indeed, the prop- 
erties of jets from active galactic nuclei can be quite dif- 
ferent on the scale of the host galaxy (see Fig. 1) and in 
the formation region (sec Fig. 2), where their transverse 
size is close to 1 pc. However, as wc mentioned above, 
the relativistic jets can be stable on all scales. On the 
other hand, in the majority of cases nonrclativistic jets 
at large distances from the star can indeed be unstable 
(see Fig. 4), so for them the treatment considered here, 
strictly speaking, can be applied only in the innermost 
parts. 

5.1 Active galactic nuclei 

Despite longstanding efforts, we know very little about 
the internal structure of jets from active galactic nuclei. 
In particular, we still do not have an answer to the key 
question of whether it is the black hole, and not the in- 
ner parts of the accretion disc, that is the central engine 
which is responsible for the black hole power [190,191]. 
In the practical sense, the main uncertainty appears in 
the determination of the particle production multiplicity 
A. Indeed, as already noted, the plasma on the mag- 
netic field lines threading the black hole horizon (which 
is needed both to screen the longitudinal electric field 
and to produce the longitudinal electric current) must be 
generated in the magnetosphere itself between two fam- 
ilies of singular surfaces. Some fraction of the plasma 
will escape the magnetosphere, while another part will 
accrete onto the black hole. Correspondingly, it is still 
unclear which quantities determine the density of mat- 
ter flowing out of the accretion disc surface (sec, for ex- 
ample, Refs [192, 193]). Moreover, it cannot be ruled 
out that the jet at large distances from the central en- 
gine will be additionally 'loaded' due to interaction with 
stellar winds from surrounding stars [194], or, for exam- 
ple, due to the 'photon breeding' effect (creation of the 
secondary electron-positron plasma by hard gamma-ray 
quanta generated in the interaction of the relativistic 
outflow with the ambient medium) [195,196]. That is 
why, the properties of jets on the scales of several kilo- 
parsecs can be significantly different from those in the 
jet formation region. 

In the region of field lines threading the black hole 
horizon several plasma generation mechanisms are cur- 
rently being discussed, in which, however, plasma is ul- 
timately generated always due to two-photon pair cre- 
ation. The one-photon conversion, which plays the lead- 
ing role in radio pulsar magnetospheres, here turns out 
to be ineffective, since the probability of pair creation in 
magnetic fields B ~ B^dd ~ 10^ G is vanishingly small. 

First and foremost, secondary plasma generation can 
be related to the direct two-photon process 7-1-7—7- 
e+ -|- e~ (see, for example, Ref. [105]), where the neces- 
sary gamma-quanta are emitted from the inner regions of 
an accretion disc. This mechanism, with a high value of 



the parameter Aagni ^ 10 -10 , was discussed in the 
pioneering paper by Blandford and Znajek [33]. How- 
ever, sufficiently high temperatures providing the neces- 
sary number of hard gamma-ray photons with energies 
above the pair creation threshold Smin = ttIqC^, and small 
free path lengths of photons are required for this mech- 
anism to be effective. Presently, the accuracy of the 
compactness parameter estimate L^agn ~ 1-100 does 
not allow one to make definitive conclusions on the effi- 
ciency of this mechanism of particle creation. 

On the other hand, such a high particle density must 
be typical for a wind outffowing from the accretion disc 
surface. Let us keep in mind that even if the energy re- 
lease related to such a wind is insignificant, it can play 
the decisive role in the matter outfiow collimation [197]. 
Here, the energy of the jet core observed at high radio 
frequencies and in gamma-rays will be associated with 
ultrarelativistic particles extracting energy from the ro- 
tating black hole. 

There is another mechanism capable of bringing parti- 
cles into the region of magnetic field lines threading the 
black hole horizon even in the absence of hard gamma- 
ray quanta. This mechanism is similar to the particle 
creation process in the outer gap of the pulsar magne- 
tosphere [198]. Indeed, the exact relativistic expression 
for the Goldreich charge density pcj takes the form [9] 

"-"-s^^'l^'S- ''''' 

In particular, near the rotational axis we simply have 



As a result, the general relativity effects cause the Gol- 
dreich density to vanish at w « fip- Therefore, a region 
quite similar to the outer gap in pulsar magnetospheres 
appears in the black hole magnetospheres. The forma- 
tion of longitudinal electric fields is also possible in this 
region, since the charge-separated plasma ffow cannot 
provide the fullfilment of the condition pe = pcj- As a 
result, it turned out that under real conditions the size 
of the acceleration region is much smaller than the sys- 
tem's size, so that the acceleration region does not affect 
the global structure of the magnetosphere [34, 35]. In 
this model, the particle production multiplicity is rather 
small: 

Aagn2 ~ 10 - 100. (123) 

Hence, we shall consider below both large and small val- 
ues of the parameter A. 

Next, it should be remembered that to explain the 
high efficiency of the energy release from the central 
engine we need to assume that the rotation parameter 
flR/c must be not too much smaller than unity. In other 
words, the light cylinder radius must not exceed the cen- 
tral engine size too much. As a result, the observed 
transverse size rjot of relativistic jets will be three-five 
orders of magnitude greater than the light cylinder ra- 
dius i?L- That is why, far from the central engine most 
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magnetic field lines must be far beyond the light cylin- 
der. 

However, according to the relationship 

^=Xr, (124) 
-Op 

where again Xr = Q.-p'Oj/c, which follows from the def- 
inition of the magnetic field components at / w /qj, 
this implies that the toroidal magnetic field will be the 
same three-five orders of magnitude stronger than the 
poloidal magnetic field. Therefore, the magnetic field 
must have a strongly pronounced spiral structure. Cor- 
respondingly, the electric field must also be three-five 
orders of magnitude larger than the poloidal magnetic 
field. At present, VLBI (Very Large Baseline Interfer- 
ometry) methods provide a lot of data on the polariza- 
tion of the innermost parts of jets [199 201]; however, 
so far it is impossible to unambiguously determine the 
magnetic field structure from the observations. 

Notice, finally, that strong twisting of magnetic field 
lines does not imply that the plasma motion will also 
occur along strongly twisted trajectories. As stressed 
above, an almost poloidal drift motion in the crossed 
electric and magnetic fields will be the principal motion 
of particles beyond the light cylinder. Thus, the toroidal 
velocity for a;^ > 1 is given by the simple relationship 

V^{Xr) = ^. (125) 

Further, if the magnetization parameter a exceeds the 
ratio rjet/Rh, the plasma Lorentz factor can again be 
estimated from the asymptotic solution (101), so that 
7 w rjci/Rh 10-^ 10^ If (7 < rjct/i?L, the effects of the 
finite mass of particles will limit their energy growth at 
large axial distances. To determine the magnetization 
parameter a, as was shown, it is necessary to know the 
particle production multiplicity A. 

As noted above, it is impossible at present to deter- 
mine the basic parameters of the outflowing plasma from 
observations. Nevertheless, some estimates can still be 
made. For example, one of the methods to determine 
the value of A is based on the assumption that the syn- 
chrotron radiation self-absorption occurs at the base of 
the jet [202]. This assumption allows one to estimate the 
particle number density [203]. On one-parsec scales, the 
characteristic particle number densities thus were found 
to be 10^-10^ cm~^, which gives 

Aagni ~ 10^° - 10^2 (126) 

If this indeed is the case, intensive secondary particle cre- 
ation near the black hole horizon should have occurred. 
Here, according to relation (56) the value of a cannot 
exceed one hundred: 

(TAGNi ~ 10^ - 10^ (127) 

which is much smaller than the ratio rjet/i?L ~ 10^ cor- 
responding to the maximum possible Lorentz factor de- 
rived from the asymptotic solution (101). Therefore, 



an almost complete transformation of the electromag- 
netic energy into particles' energy must occur in the pro- 
cess of collimation. Notice that in this case, although 
7max ~ (T ^ 10^-10"^ exceeds the particle energies that 
are required to explain the apparent superluminal mo- 
tion effect, it is still insignificant. Here, almost all the 
energy flux in the jet will be related to the flux of acceler- 
ated particles. For a ^ 10^ lO'', the radius rp ^ a^^^R^ 
(66) of the fast magnetosonic surface must be smaller 
than one hundred radii from the central engine, which is 
10^*-10^^ cm. Thus, the flow in the jet must indeed be 
supersonic. 

Now making use of expressions (89), (95), and (96), we 
arrive at the conclusion that in this case all critical mag- 
netic fields for the reasonable value of 7i„ ~ 10 are larger 
than Bcxi ^ 10^^' G corresponding to the ambient pres- 
sure (for convenience, all parameters discussed here are 
listed below in Table 2). But this means that a denser 
core must exist in the center of the jet, and the subsonic 
flow near the axis will not be formed. At last, the ejec- 
tion rate of electron-positron pairs, N = ttRq ■ Ancjc, 
can be estimated, which, as can be easily checked, is 
determined using the simple relation 

N^x(^^y\ (128) 

As a result, we have iVAGNi ~ 10^^ particles s^^ (here- 
inafter we set M = lO^M©). In other words, about 
10^^ electron-positron pairs will be injected into host 
galaxy over the time of the active life of a galactic nu- 
cleus, T ^ lO*" years. This number, incidentally, is 
quite sufficient to explain the intensity of the annihi- 
lation line emitted from the Galactic center, which, as is 
well known, requires about 10^^ annihilations per second 
[204]. 

If the multiplicity factor of the secondary particle cre- 
ation is small, AagN2 ^ 10 100, the inner structure of 
the jet must be significantly different, since now all crit- 
ical fields are below the value of 10~^ G corresponding 
to the ambient pressure. Here, one obtains 

CTAGN. ~ 10^° - 10^^ (129) 

so that the plasma Lorentz factor, according to the 
asymptotic solution (101), can be as high as approx- 
imately 10^ 10"'', and a subsonic fiow region must be 
formed in the center of the jet. In this case, the energy 
flux is Poynting-dominated. Correspondingly, the dense 
core will not be formed, so that both the poloidal mag- 
netic field and the plasma density are weakly dependent 
on the axial distance. The electron-positron pair injec- 
tion rate will be much smaller: iVAGNa ~ 10^^ particles 
s~^. However, in this case, too, the fast magnetosonic 
surface radius rp ~ 10^''-10^^ cm will be smaller than 
the jet transverse size. Notice that here there is no di- 
rect contradiction with observations, since, as has been 
noted, the drift motion of particles will be directed al- 
most along the poloidal magnetic field. This particle 
motion does not produce synchrotron radiation. Hence, 
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one should be cautious in using the standard synchrotron 
radiation formulas to estimate the value of the magnetic 
field and the lifetime of relativistic particles. 

Wc should keep in mind that, when estimating ra- 
diation from relativistic jets, one usually assumes that 
an approximate cquipartition between the energy den- 
sities of particles and the magnetic field takes place 

(-^(com)2/gj^ _ ^(com)^(com)^^2-) ^j^^ ^.^^^ ^^.^^-^^ 

plasma. The parameters we introduced allow us to ob- 
tain simple relationships for all quantities in this refer- 
ence frame. In particular, it is easy to show that the 
characteristic Lorentz factor of particles in the plasma 
rest frame is simply the following: 

7 

where 7 is the bulk Lorentz factor of the hydrodynamic 
flow. On the other hand, one finds S^™'") w {xr/j)Bp. 
Consequently, we have i?^'^"™) « Bp in the effective ac- 
celeration region (7 = Xr), while for the asymptotic so- 
lution (102) we obtain > Bp. 

5.2 Micro quasars 

If the operation of the central engine in microquasars 
indeed can be described by the MHD model considered 
here, it is again possible to assume that the observed sub- 
relativistic jet velocities must correspond to a on the or- 
der of 3-10. Then, according to relation (56), we should 
conclude that the particle production multiplicity in the 
microquasar magnetosphere must be fairly large: 

A^QSO ~ 10'°. (131) 

Such a high value is also supported by the compact- 
ness parameter Za./^QSO ~ 10"*. Then, the electron- 
positron pair ejection rate should be iV^qso ~ 10*'^ par- 
ticles s^^. Finally, the large values of the critical fields 
(Smin ~ 10"^ G, Beq ^ 10^ G, B„ ~ 10*^ G) indicate 
that a dense core must exist in the jet center, and the 
subsonic flow near the rotational axis docs not form. 

On the other hand, if a is indeed not too high, the fast 
magnetosonic surface (rp ~ cr^^^Rh) must lie close to the 
light cylinder, i.e., at distances of order 10^-10® cm. This 
scale is much smaller than the distance from the central 
engine to the companion star in binaries involving a mi- 
croquasar. Hence, we can conclude that the interaction 
of the stellar wind and matter ejected from the micro- 
quasar occurs in the supersonic regime. 

5.3 Sources of cosmological gamma-ray 
bursts 

Now let us see which parameters can characterize jets 
outflowing cosmological gamma-ray bursts. It should 
be recalled that, according to one of the most popu- 
lar models, a rapidly rotating central engine (magnetar, 
black hole) is immersed in the progenitor massive star 



[40]. That is why, the ambient pressure for the jet in 
its formation region is provided not by the surround- 
ing medium with very small pressure, but by the mat- 
ter of the massive star itself (the equivalent magnetic 
field -Bext ~ 10^-10^ G). Notice also that in cosmological 
gamma-ray bursts there can be one more mechanism of 
the electron-positron pair creation, related to neutrino 
annihilation. Such neutrinos can be copiously created 
during supercritical accretion onto the collapsing stellar 
core [205]. 

The starting point that can shed light on the physical 
conditions inside the central engine can be the character- 
istic particle Lorentz factor 7 ~ 300, which in this case 
can be naturally related to the magnetization parameter 
a: 

(7GRB ~ 10^ - 10^. (132) 

If the condition 7 <C o" is satisfied, the total energy re- 
lease from the central engine would be unrealistically 
high. Now, using Eqn (56) to estimate the plasma gen- 
eration multiplicity, we obtain 

Agrb ~ lO'^ - 10^^ (133) 

Such a huge value unambiguously evidences that the 
particle creation efficiency must be high enough. In- 
deed, formula (57) gives very large compactness parame- 
ter la ~ 10^^. Correspondingly, for the electron-positron 
pair ejection rate we find A^grb ~ 10^^ particles s~^. 

Next, the very small sizes of the central engine to- 
gether with the moderate value of the magnetization 
parameter a shows that the fast magnetosonic surface 
radius rp (66) should not exceed 10''-10^ cm, which is 
significantly smaller than the size of the progenitor star. 
Consequently, the matter outflow in the jet becomes su- 
personic before it exits the star. Finally, expressions (89) 
and (95) for the characteristic magnetic fields -Bmin and 
Beq indicate that they are in the range of 10^-10^ G, i.e., 
their pressure is comparable with that inside the progen- 
itor star. Therefore, the jet transverse size will indeed 
be sufficient to accelerate particles up to energies -f ^ a. 

Notice at last that the condition 7'!9 ^ 1 is certainly 
not satisfied for gamma-ray bursts, since in that case the 
jet spread angles would only be 0.1°, while observations 
indicate that -d ^ 1-10° [206]. Such a flow can also be re- 
alized. For example, it was shown in paper [207] (see also 
Ref. [208]) that in the model of an infinitely extended 
progenitor star, where the ambient pressure decreases 
graduately according to a power law, the acceleration 
turns out to be not very effective in comparison with a 
more realistic model in which the ambient pressure be- 
yond the star is assumed to be low. As it has turned 
out, in both cases the flow corresponds to a weakly col- 
limated flux with 1 < fc < 2, where the particle energy 
follows the asymptotic behaviour 7 « [Rc/zuY^'^ (102). 
In the former case, however, the radius of curvature Rc of 
the magnetic surfaces, which is determined by the pres- 
sure decrease law inside the progenitor star, turns out to 
be sufficiently small, which precludes plasma from being 
effectively accelerated. Beyond the star, magnetic field 
lines straighten up (and hence the radius of curvature 
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increases), which leads to a more effective acceleration. 
As we see, the simple analytical asymptotic solutions ob- 
tained above allow easy interpretation of the numerical 
experiment. 

5.4 Radio pulsars 

Radio pulsars, undoubtedly, only indirectly relate to the 
topic considered here, since their magnctosphcrcs are 
certainly not axisymmetric and stationary. It is not then 
surprising that jets, as we have noted, are observed only 
from two energetic pulsars. Nevertheless, many points, 
which were possible to clarify in the theory of the neu- 
tron star magnctosphcrc, undoubtedly allow us to shed 
light on the nature of other compact objects, too. 

First and foremost, it should be noted that basic pa- 
rameters characterizing the pulsar wind are well known 
inside the light cylinder. This is due to our good knowl- 
edge of the process of plasma creation near the magnetic 
polar caps of a neutron star. Numerous calculations 
have shown that the general properties of the secondary 
electron-positron plasma flowing out from the magne- 
tosphere turned out to be only a little sensitive to the 
details of the acceleration region structure. For most 
models [178, 180, 209, 210], both the density and the 
energy spcctriim of the plasma flow ing out arc universal 
enough. Therefore, with certainty we can say that the 
plasma flowing along open field lines in the pulsar mag- 
netosphere includes both the primary particle beam with 
energy £ « lO'' MeV and density close to the Goldre- 
ich density and the secondary electron-positron compo- 
nent. Its energy spectrum with a good accuracy exhibits 
a power law dependence 



£ oc £ 



(134) 



ranging from ^min ~ 10-100 MeV to 5max ~ 10^ MeV. 
The total density of the secondary plasma for ordinary 
pulsars exceeds the Goldreich density by lO'^-lO'' times. 
And only for the most energetic pulsars can the multi- 
plicity factor reach 10^. 

Thus, the parameter A for radio pulsars is determined 
quite reliably: 

ApsR ~ 10^ - lO''. (135) 



Now, making use again of relation (56) we obtain 
(TPSR - 10' - 10^ 



(136) 



And only for the most energetic pulsars do we find 
crpsR ~ 10^-10^. Thus, the condition 7in » cr^/^ is sat- 
isfied in the vast majority of pulsars, which corresponds 
to a slow rotation of the central object [99]. Indeed, 
the equality 7in = cr^^^ can be written out in the form 
P = Per, where 
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For fast rotation {P 4C -Per)) the particle energy signif- 
icantly increases when particles approach the fast mag- 
netosonic surface, whereas for slowly rotating objects 
(P S> Per the particle energy remains practically un- 
changed. Their further fate, as we have shown above, 
depends on whether the flow intersects the light surface 
or not. 

5.5 Young stars 

In conclusion, we discuss the main parameters which 
characterize nonrelativistic jets from young stars. It 
should be recalled that in this case the nonrelativistic 

magnetization parameter an an ~ (rjp/f^cr)^ (53) plays 
the key role. Under the condition CTn ^ 1 (^^f >• f^cr), 
the electromagnetic energy flux near the central engine 
will be much greater than the particle energy flux; be- 
yond singular surfaces, as stressed above, particles must 
carry a significant fraction of the total energy. 
As the critical period Per = 27r/r2er7 namely 
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(fd is the inner radius of the accretion disc) is close to 
the spin periods of young stars, in the region of magnetic 
field lines coming out of the surface of the star we have 
(Tn ^ 1. On the other hand, the period Per is one two 
orders of magnitude larger than the rotation periods in 
the inner regions of accretion discs, Pb = 27r{GM/r^y^'^ , 
so that for the corresponding magnetic field lines one 
finds (Jn ~ 10-1000. That is why, the inner parts of the 
accretion disc, not the central star, must play the role 
of the central engine rotor. As noted above, there are 
observational evidences of this being the case [76] . 
Next, from relation (66) we obtain 



Vin 1/3 



(139) 



so that rp ~ 10-30 P. Thus, the distance to singular 
surfaces exceeds the size of the star by 10-30 times, 
but it is the same 10-30 times smaller than the trans- 
verse size of jets. Consequently, the flow inside the jets 
must be supersonic and, hence, the longitudinal current 
for these objects again must be derived from the crit- 
ical conditions on the singular surfaces. Correspond- 
ingly, the radius of the jet core must be on the order 
of 7'coro = Vin/^ ^ 0.1 a.u., the jet magnetic field should 

be as high as Peore ~ 0.1(rircore/t'm)^Pm ~ 10"^ G, 

and the particle number density must range from 10^ to 
10^ cm-3. 

Finally, we note that the existence of the integrals of 
motion allows us to obtain direct information about the 
plasma outfiow region. For example, if the radial and 
longitudinal velocities of the flow in the jet are known at 
the axial distance r± (and such observations, as noted 
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above, have already been performed for several young 
stars [72, 73]), it is possible to estimate the distance Vst 
from the central star at which the corresponding force 
line is 'anchored' in the accretion disc [211, 212]: 
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It is obvious that this scale is much larger than the size 

of the star, so at present it has been possible to resolve 
only outermost regions of the outflow. 

6 Conclusions 

Thus, significant progress has indeed been recently 
achieved in the understanding of the natiue of jets ob- 
served in different classes of astrophysical sources. This 
became possible because the analytical approach allowed 
sufficiently simple relationships between physical param- 
eters characterizing the outflows to be found, and the nu- 
merical modeling (in which, we should keep in mind, the 
setup of the problem itself has been significantly different 
from the analysis of stationary equations) has confirmed 
these relationships. 

The most important result of the analytical theory in- 
cludes the understanding of the role of key dimensionless 
parameters. For clearness, they arc listed in Table 2 (for 
active galactic nuclei we set M = 10^ Mq). It turned out 
that the knowledge of these parameters allows us to esti- 
mate many characteristics of jets, including the fraction 
of energy carried by particles, the plasma Lorentz factor, 
the electron positron pair injection rate N, and the com- 
pactness parameter la, as well as to determine the main 
parameters of the internal structure of jets. The determi- 
nation of these parameters from observations would be a 
significant breakthrough in our understanding of physi- 
cal processes which are underway in active astrophysical 
sources. 

Next, we have shown that many properties of rela- 
tivistic and nonrelativistic jets are significantly different 
from each other. For convenience, we also collect them 
together in Table 3. As can be seen, one should be very 
cautious when applying the results, which were obtained 
for nonrelativistic jets, to ultrarelativistic flows. More- 
over, the asymptotic relations formulated above clarify 
many results obtained by means of numerical simula- 
tions. 

The limited space of the present review did not allow 
us to discuss in detail many important issues. In particu- 
lar, we only briefly discussed the stability of jets. Finally, 
here we have no space at all to discuss the nonstation- 
ary performance of the central engine (most papers have 
recently started focusing exactly on this topic) or the 
proper radiation of jets. Nevertheless, we would like to 
hope that questions addressed in this review will be use- 
ful for future studies of relativistic and nonrelativistic 
outflows observed in many astrophysical objects. 



Table 2: Parameters of jets outflowing from relativistic 
compact objects. All values are given to an order of 
magnitude. 





AGNi 


AGN2 


mQSO 


GRB 


PSR 


a 


100 


1012 


10 


10=^ 


10^ 


A 


1012 


100 


lOlo 


1014 


10^ 


7in 


10 


10 


10 


10 


100 


la 


1-100 


1 100 


104 


lOi^ 


10-'^ 


Wpart/VFtot 


1 


10-9 


10-5 


1 


10-2 


7 


10-100 


lo-i-io^ 


10^ 


300 


10^ 


N, s-i 


1049 


1039 


1033 


1053 


1032 


-Bmin, G 


10-2 


10-12 


104 


10^ 


fo-« 


Bcq, G 


fo-4 


10-24 


10^ 


10^ 


10-8 


BcT, G 


10-1 


10-11 


lO'^ 


101° 


10-2 
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7 Appendix 

In the Appendix, we give for reference the complete set 
of equations of the Grad-Shafranov method written out 
in themost general case, i.e., for axisymmetric stationary 
flows in the vicinity of a rotating black hole. First of all, 
we keep in mind the basic relations for the Kerr metric of 
a rotating black hole. In the Boyer-Lindquist coordinates 
t, r, 6, and ip, it assumes the form 

ds2 = -a'^dt^ + Qikidx' + I3'dt){dx'' + jS'^dt), (A.l) 

where the quantity 

a=^VA (A2) 



is the gravitational redshift, and the vector ^ is toroidal: 
= = 0, = -cj. (A3) 

Here 

9n Mr 

(A.4) 



2aMr 



is the so-called Lense-Thirring angular velocity. Finally, 
M and a are the mass and specific angular momentum 
of the black hole; (a = J/M), respectively. In addition, 
we introduced the standard notations 
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Table 3: Main differences between rclativistic and nonrelativistic jets. 



Relativistic fiow 


Nonrelativistic fiow 


Longitudinal current is close to the Goldreich current 


Longitudinal current is much larger than the 

Goldreich current 


For strongly magnetized flow (cr ^ 1), fast 
magnetosonic surface near equatorial plane is located 
cr^/^ times farther than Alfvenic surface 


Fast magnetosonic surface is located near Alfvenic 
surface 


On fast magnetosonic surface, the particle energy flux 

is much smaller than the electromagnetic energy flux 


On fast magnetosonic surface, the particle energy flux 

is close to the electromagnetic energy flux 


Beyond Alfvenic surface, the electric field is close 
in magnitude to the magnetic field 


Beyond Alfvenic surface, the electric field is much 
smaller than the magnetic field 


Proper coUimation is impossible for both strongly and 
weakly magnetized outflows 


Proper collimation becomes efl'ective for strongly 
magnetized outflows 


The dense core in jets can be formed only under 
sufficiently small ambient pressures. Magnetic flux 
in the core is much smaller than the total flux 


The dense core is always formed in jets. Magnetic flux 
in the core is a significant fraction of the total 
magnetic flux 


Cylindrical flow with subsonic core is possible 


Cylindrical flow with subsonic core is impossible 



Here, in all relativistic expressions we use the units in 
which c = G = 1. Finally, it is important that the three- 
dimensional metric gtk in formula (^.1) be diagonal: 



9rr — ^ , 
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(A.6) 



As for the flat space limit, it can be easily obtained by 
passing to the limit a — >• 1 and w ^ 0. 

As is well known, for calculations it is convenient to 
introduce a special reference frame, the so-called ZAMO 
(Zero Angular Momentum Observers) [131], which has 
the following properties: 

• ZAMO observers are located at constant radius 
r = const, 6 = const but rotate with the Lense 
Thirring angular velocity dip/dt = lo; 

• for ZAMO, the four-dimensional metric gaB is diag- 
onal, with its three-dimensional part gik coinciding 
with Eqn (A.6). 

Below, all vectors will be written out in this reference 
frame. In particular, the operator Vj, means the covari- 
ant derivative in the three-dimensional metric (A.Q). 

As a result, the electric and magnetic fields can be 
conveniently written as 



i> = 



21 



E : 



2ttzu aw 

Op — w „ T 
-V'l', 



6(3, 



{A.7) 



2Tra 

respectively, and the four-velocity of matter is written 
as 

u = — B -I- 7(s2f - w)— e<3, 
an a 



{A.9) 



where 7 = — v"^ is the Lorentz factor of matter, and 

the subscripts with a cap over them correspond to physi- 
cal components of vectors. The quantity fip remains the 



integral of motion. In turn, integrals of motions E and 
L will be written out as 



E = £;(*) 



~2^ 



+ l^r]{aj + ojzuu^), (A.IO) 



L = L{^) = — + lirjwu^. {AM) 

ZTT 

Next, the relativistic Bernoulli equation 7^ — u| = Up-|-1 
takes the form 

^ = 64;?^^+"W, (A12) 

where the Alfvcn factor is 



(A.13) 



and 



K = a'^w^{E - ClpLf [a2 _ (Qp _ u))^^^ _ 2M^' 



(A.14) 



This equation defines the Alfvenic Mach number M, 
where 

2 47r?yV 



(^.15) 



Now, making use of relations (A.12)-(A.14), which can 
be recast in the form (V*)^ = F{J\A'^ ,E,L,ri,Q.F, jj), 
where 



647r4 K 647r4 2222 
F = —r-rr— r-rra ZU T] , 



we obtain 



(A.16) 
(A.17) 
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where 



Us 



1 /4 V' F 



1 {E - QfL) (Of - (^)ru'^ - LM^ 
zufXT] — (rip - 'jjf-w'^ — M? 



(y4.26) 



(V*) 



2 (V*)2' 



(A18) 



Here, the operator V'j, acts on all quantities but . 
The quantity Vjl./i must be determined from the relation 
[97] 



1 (dP 



n \ ds 



+ T 



VfeS, 



, 2cg VfcT? , 1 

(A.19) 

where = 1/ ^{dP/dn)s is the speed of sound, and 
s is the entropy. In turn, the denominator D can be 
rewritten in the form 



D 



Bl 



1 A 



{A.20) 



As for the Grad-Shafranov equation, in the compact 
form it can be written out as [96, 97] 



Equations {AA2) and (A.24)-(A.26) represent alge- 
braic bounds which allow the determination, albeit in 
an indirect form, of all characteristics of the flow from 
the given poloidal field Bp (i.e., from the known po- 
tential 5') and five integrals of motion. It should be 
emphasized that for a nonzero temperature they are ex- 
tremely lengthy, mainly due to the need to resolve equa- 
tion (A. 19). In the case of cold plasma (s = 0, i.e., 
fj, = const), Bernoulli equation (A. 12) becomes a fourth- 
order algebraic equation with respect to M^. As shown 
above, this fact often allows analytical asymptotics to be 
found. 

In the cylindrical case, the second-order Grad- 
Shafranov equation can be conveniently reduced to the 
system of two ordinary differential equations of the first 
order for the magnetic flux 'i!{w) and the Mach number 
A4^. The equation for the Mach number has therewith 
the form [146] 



a 



1 



aw 



+ 



, dfip 
"d^ 



(^.21) 
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dr Ar^^'^rfc'^ 
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{e'f 
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d* de' 



yp'rfd^ dr d^* 



647r4 1 d (G\ 1 dr? 3 ^ ds „ 

a^vj^ 2M^ d-^ \A J r/d* d* 

where 

G = a^w'^{E - n^Lf + a^JveL^ - M^zu^{E - ujLf. 

{A.22) 

Now, expanding terms VfcA^^ in Eqn (A. 21) according 
to definitions {AA7)-{AA^), we finally arrive at 



+M- 



d* dOp 



dr 

A 
n 



-M' 1 



2A- 



cl \ I dr] 
dr rj d^" 



dP\ 



M^d^ds_ 
II dr d*' 



(A.27) 



where e' = E — VIyL. The equation for the magnetic flux 
will coincide with Bernoulli equation (A. 14). Finally, 
the force-free pulsar equation takes on the form 



+ 



a2n72(V*)2 
A 1 



o o k 
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vo dw 



a2n72(V*)2 2D 
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{A.28) 



Idr; ,„Q^ds 
-167r^un-— f - len^nT— = 0, 
r?d* d* 



{A.23) 



wlicire is the Laplace operator. Its generalization 
to the force-free black hole magnetosphere is written as 
[131] 



where again the gradient acts on all quantities but 

Ai'^, and the derivative d/d'^ acts only onh the inte- 
grals of motion. Formula (j427) determines in the most 
general form the equilibrium equation for the magnetic 
surfaces. Finally algebraic relations have the form 



1^ f a 



1 



(Qp - co)^w^ 



+ - 



fip — U) 



167r2 dl ^ 
+ =0. 



{A.29) 
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(yl.24) 
(A.25) 



These equations are elliptical in all the space, and so they 
require boundary conditions on the integration region 
boundary or on the black hole horizon. 
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